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High-Speed Testing in the Southern 


SUMMARY 


The basic performance and flow characteristics of this variable 
density, high-subsonic-speed wind tunnel are described. The 
problems of model support and of special mounting systems for 
wind-tunnel tests at Mach Numbers close to 1 are discussed in 
detail. In particular, the Lockheed “Bump” method of testing 
through the transonic region in a wind tunnel is considered. 
Calibration, operatiop, and data reduction procedures required 
for obtaining and interpreting experimental data at Mach Num- 
bers near 1 are discussed, especially the effects of compressibility 
on blocking and other tunnel-wall corrections. Pertinent ex- 
perimental results from a number of recent tests are presented. 


INTRODUCTION 


# See MAIN CONSTRUCTIONAL FEATURES of the South- 

ern California Cooperative Wind Tunnel (CWT) 

* have been described'!; the aerodynamic design was dis- 
cussed at the August meeting of the Institute of the 
Aeronautical Sciences in 1945. Since the latter presen- 
tation was never published, certain of the more signifi- 
cant elements will be referred to again in the present 
paper, whose main objective is the discussion of some 
aspects of the problem of high-subsonic- and transonic- 
speed wind-tunnel testing. 

The CWT has now been in essentially routine opera- 
tion for about a year and the major portion of the tests 
carried out have been concerned with the high-subsonic- 
speed range. Not only are questions of model support 
and wind-tunnel wall interference extremely difficult 
and complex in this region, but even the fundamentals 


Presented at the Annual Summer Meeting, I.A.S., Los Angeles, 
August 7-8, 1947. 

* The authors wish to express their appreciation of the assist- 
ance of many members of the staff, especially that of Keirn Zebb 
and Donald Francis, who made important contributions to the 
paper. 

T Director, Assistant Director, and Coordinating Supervisor, 
respectively. 

t Owned jointly by the Consolidated-Vultee, Douglas, Lock- 
heed, and North American aircraft companies, and operated by 
the California Institute of Technology. 
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of the simplest transonic flows are not understood. 
Indeed, the field of transonic aerodynamics furnishes 
some of the most difficult of the unsolved problems of 
modern fluid mechanics. Accordingly, there has been, 
and continues to be, an extreme need for systematic re- 
search investigations designed to clarify the basic 
phenomena involved. 

Unfortunately, the pressure for tests of specific 
models has been so great that such systematic investi- 
gations have not as yet been possible at the CWT. 
In many cases even essential calibrations have had to 
be fitted in between other tests and have often had to 
be dropped before the desired scope or accuracy had 
been attained. Asa consequence, the present paper will 
inevitably appear incomplete and perhaps unsystem- 
atic. In this connection the authors can only state their 
keen awareness of the deficiencies and gaps and express 
the hope that, in view of the current interest in the 
problem, their experience may be of some assistance to 
other workers in the field. 


TUNNEL DESCRIPTION 


Fig. 1 gives a general view of the more important ele- 
ments of the tunnel, which is normally run with the 
large circular door closed and sealed. Under these 
conditions it can be operated at any internal pressure 
in the range between one-tenth and three atmospheres. 
The working or test section is a rectangular box 18 ft. 
long, 12 ft. wide, and 8'/2 ft. high, with large corner 
fillets. Gates at each end of the test section can be 
closed so as to seal off the working area from the rest of 
the tunnel and permit the former to be brought quickly 
to atmospheric pressure for model access without dis- 
turbing the pressure in the major portion of the tunnel 
circuit. It should be mentioned that the hydraulically 
operated and pneumatically sealed gates and main door 
operate so rapidly that it is practically as simple to run 
the tunnel at any nonatmospheric pressure in the design 
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Fic. 1. Cutaway drawing of the CWT. 


range as at atmospheric pressure. It follows that tests 
can be made at whatever pressure seems most desirable 
without any appreciable sacrifice in operating efficiency. 
The flexibility thus achieved has proved to be one of 
the tunnel’s most important characteristics. 

Air is circulated through the tunnel by a 32-bladed, 
two-stage propeller-type fan driven by electric motors 
with a 12,000 hp. maximum output. A coupling be- 
tween the two stages permits the downstream one to be 
disconnected from the drive shaft for high-pressure 
operation, and the blades of the two stages are sepa- 
rately controllable with the fan stationary or rotating. 
This flexibility in control, combined with variable- 
speed drive motors, makes it possible to operate over 
the full tunnel pressure range with satisfactory fan ef- 
ficiency. 

A radiator at the corner vanes upstream of the work- 
ing section removes enough heat so that, under condi- 
tions of protracted full-power operation, the stagnation 
temperature remains below 125°F. Under usual oper- 
ating conditions this temperature seldom exceeds 
100°F. 

One of the most important features of the tunnel, 
which is of particular significance in connection with 
high Mach Number tests, is the fact that there are 
three interchangeable working sections, any one of 


which can readily be rolled into the tunnel so that the 


model support system is rigidly attached to the metrical 
or balance system. Two of these ‘‘carts,’’ as they are 
called, are designed for strut-type model supports 
(utilizing one, two, or three struts), while the third in- 
cludes a ring frame surrounding the working section so 
that wing models spanning the throat can be mounted 
from flush face plates in either a horizontal or vertical 
plane (12 and 8'/2-ft. span, respectively). This pro- 
vision makes possible the use of-a wide variety of dif- 
ferent types of model support systems and has proved 
of special value for investigations close to M = 1.0. 

The metrical elements of the balance system are 
Tate-Emery hydraulic capsules connected to the bal- 
ance indicators by high-pressure oil lines. This hy- 
draulic system has rapid response and intrinsically high 
damping. Experience has shown that it is especially 
suitable for high Mach Number tests where flow in- 
stability and model vibration lead to forces that are 
fluctuating with extremely large amplitudes and fairly 
high frequencies. Even close to choking conditions 
where the entire tunnel structure, including the balance 
system attached to it, is vibrating severely, the balance 
readings are still steady and apparently give satisfac- 
tory mean values of the fluctuating forces. 

The most serious difficulty with the balance system 
arises from the fact that it had to be designed to measure 
maximum lift and drag forces of 30,000 and 5,000 Ibs., 
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respectively, while experience has demonstrated that 
for measurements close to M = 1 it is essential that the 
models be small and the tunnel pressure be well below 
atmospheric. The result is that the forces to be meas- 
ured are small and the accuracy possible with the mas- 
sive metrical system becomes marginal or inadequate. 
The following numerical values correspond to a typical 
high Mach Number investigation carried out on an air- 
plane model with 1-sq.ft. wing area at a tunnel pressure 
of one-half atmosphere. These conditions make pos- 
sible the attainment of a Mach Number close to 1 
through the use of somewhat less than full available 
power. The dynamic pressure is about 450 lbs. per 
sq.ft., and the uncertainties in the three major coef- 
ficients are: 


AC, = 0.01, AC,, = 0.003, AC, 0.0010 


If the pressure is held constant, these uncertainties in- 
crease by a factor of about 5 at M = 0.35. On the 
other hand, if the pressure is increased and the tunnel 
power is held roughly constant as M decreases, the 
errors decrease and reach values about one-third of 
those listed at M = 0.35. 

It appears that a 1-sq.ft. model is about the smallest 
for which the balance system is at all suitable. Ac- 
cordingly, special strain gage balances have been de- 
veloped for use with certain small models. In some 
cases it is possible to install such a balance inside the 
model, thus eliminating the tares (but not necessarily 
interferences) arising from the model supports. It is 
expected that this strain gage technique will be inten- 
sively developed during the next few years. 

One difficulty that has plagued all high-speed wind 
tunnels is ‘‘sandblasting’’ of models by solid particles 
carried by the air stream. This problem has been at- 
tacked at the CWT by mounting along the vanes in the 
corner downstream of the fans a large mesh wire screen 
covered with fine cheesecloth. There is usually a cer- 
tain amount of oil vapor in the tunnel coming from the 
fan and balance systems, and this makes the cheesecloth 
an excellent dust catcher. Following a model installa- 
tion or any considerable work done in the tunnel cir- 
cuit, some dust and solid particles remain in spite of 
the most careful cleaning. However, after a couple of 
runs at high speed, practically all of the solid particles 
have been collected by the cheesecloth and no further 
sandblasting of the model occurs at even the highest 
Mach Number. 


Basic FLow CHARACTERISTICS 


The true, dyed-in-the-wool wind-tunnel operator is 
always a little unhappy at the insistent demands that 
are made upon him for actual model tests, since he is 
never completely satisfied with the calibrations that he 
has been able to complete on his enormously complex 
instrument. Adequate calibration of the CWT has 
been especially difficult because of the complications in- 
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Fic. 2. Mach Number distribution at two radii in horizontal 
plane at center of working section. 


herent in the operating flexibility designed into the 
tunnel, because of its primary function of assisting in 
the solution of the aircraft industry’s ever present and 
urgent immediate problems, and, perhaps most imp or- 
tant, because of the intrinsic difficulties in the transonic- 
speed régime with which it has been most concerned. 
The more important flow characteristics that have been 
determined during such calibrations as have been pos- 
sible are presented in this section. Each of the investi- 
gations for which results are reported has suggested 
many further studies that should be made but which 
could not yet be pursued. It is hoped that it may be 
possible in the future to carry at least some of them 
much further. 

The velocity distribution in the region near the center 
of the working section, normally occupied by a model, 
was determined by surveys with a high-speed pitot- 
static tube mounted on a thin, swept-strut support. 
The survey was conducted by setting the pitot tube at a 
fixed radius from the center point of the working section 
and rotating it in a horizontal plane around a complete 
circle by means of the model yawing mechanism. The 
pitot tube and strut were kept accurately aligned with 
the tunnel axis throughout the traverses by a counter- 
rotating mechanism built into the model suspension 
system. Typical results for three nominal Mach Num- 
bers and two radii are shown in Fig. 2. It appears that 
the velocity distribution over the central portion of the 
working section is satisfactorily uniform. 

The nominal or “index” Mach Number, M,, used 
throughout the paper is the Mach Number at the center 
of the working section, with only the measuring instru- 
ment installed, for a given ‘‘venturi pressure difference,” 
Ap,. This pressure difference is indicated by two piezom- 
eter rings, one located in the large-diameter entrance 
section and the other a few feet upstream of the work- 
ing-section leading edge. 

The basic calibration or determination of M, was 
accomplished by means of the “‘static pipe.’’ This is a 
polished steel, circular tube approximately 5 in. in 
diameter and 45 ft. long, containing 34 pairs of care- 
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Fic. 38. Axial Mach Number distributions from static pipe and ceiling orifices 


fully smoothed static pressure orifices spaced at 1-ft. 
intervals. The tube extended upstream into the slow- 
speed entrance section and downstream into the dif- 
fuser. All supports were upstream or downstream of 
the high-velocity working-section region (Fig. 3). 
Pressure orifices were also installed along axial lines in 
the ceiling of the working section and diffuser. The line 
of orifices was located approximately 2 ft. away from 
the ceiling centerline in the working section and close 
to the centerline in the diffuser. 

Axial Mach Number distributions computed from 
the observed pressures at these sets of orifices are plot- 


ted in Fig. 3. Unfortunately, the two sets of data were 
not taken simultaneously, so that the index Mach 
Numbers for pipe and ceiling observations are not al- 
ways identical. The Mach Number is uniform over 
the downstream portion of the working section and, 
although choking occurs at the beginning of the dif- 
fuser, the Mach Number at the center of the working 
section is only slightly less than one at the choking 
condition. This indicates that the slight downstream 
increase in working-section area, which was designed to 
compensate for boundary-layer growth, was very 
nearly correctly chosen. The sharp increase in ceiling 
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Fic. 4. 


velocity due to local flow curvature at the junction be- 
tween working section and straight-walled diffuser is 
clearly shown. 

In going upstream from the working-section center 
the velocity decreases appreciably, the effect becoming 
large at high Mach Numbers. This indicates that the 
downstream boundary-layer growth was not ade- 
quately compensated in this region. This effect is 
shown more clearly in Fig. 4 where is plotted axial dis- 
tribution of pressure, p,, referred to static and dynamic 
pressure at the working-section center (p,, g.). The 
area distribution that would lead to the pressure dis- 
tribution, assuming one-dimensional adiabatic channel 
flow, has been computed for each of the experimental 
curves. The results are also plotted in Fig. 4 in the 
form of a dimensionless difference between area at an 
arbitrary axial position, A,, and area at the channel 
center, A,. Except for the region near the beginning 
of the diffuser, all of the area curves coincide with satis- 
factory accuracy. It appears that the effective area at 
the beginning of the working section is a little less than 
| per cent larger than that for the downstream portion. 


Axial pressure and effective area distributions from static pipe measurements. 
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Pressure coefficient distribution along diffuser, theoreti- 
cal and experimental. 


Fic. 5. 


It is planned to decrease the upstream area by building 
up the corner fillets in accorda:ce with the area distri- 
bution curve so as to give an essentially constant static 
pressure throughout the entire working-section length. 
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Fic. 6. Boundary-layer profiles at upstream and downstream 
positions in the working section. 


The diffuser data of Fig. 3 have been replotted in Fig. 
5 so as to give a simple visual indication of diffuser 
efficiency. The theoretical pressure coefficient, C,,, 
computed from the axial area distribution on the basis 
of one-dimensional frictionless flow of an incompressible 
fluid, is plotted as the solid curve. Each experimental 
pressure coefficient has been multiplied by the ratio of 
the theoretical coefficient at M = 0 to that at the Mach 
Number of the test. If there were no diffuser losses 
all the experimental points, as so modified, would fall 
on the theoretical 4 = O curve. It appears that the 
diffuser efficiency is extremely high and essentially in- 
dependent of Mach Number, except for the highest 
test Mach Number where a shock wave had already ap- 
peared at the entrance to the diffuser. 

Measurements of boundary-layer profiles have been 
made at various points on the working-section ceiling 
at a number of Mach Numbers, using total pressure 
combs. Two sample sets of data are presented in Fig. 
6 for upstream and downstream locations in the working 
section. It appears that the boundary-layer thickness 
increases from about 2 in. to about 6 in. in going down- 
stream and that it varies relatively little with Mach 
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‘Fic. 7. Scale and intensity of turbulence from hot-wire meas- 
urements. 


Number. Another set of measurements (not reported 
here) was taken at a point on the floor directly below 
the center of the working section. The results indicate 
a boundary-layer thickness of approximately 3 in., 
essentially independent of Mach Number, at this section 
where models are usually located. These values of 
boundary-layer thickness are slightly larger than had 
been anticipated on the basis of smooth-walled model 
tests but are not considered excessive in view of the 
tee slots, plate joints, and other irregularities on the 
working-section surfaces. 

Turbulence measurements have recently been made 
using a tungsten hot-wire and special equipment de- 
veloped by Dr. Liepmann’s group at the GALCIT. 
The hot-wire was located at the center of the working 
section, and observations were taken at one half and 
one atmospheres tunnel pressure. The tungsten hot- 
wire made it possible to extend certain of the readings 
up to M = 0.7. The results, both as to intensity and 
scale? are presented in Fig. 7. The lateral scale, L, 
was computed from observations of the longitudinal 
scale obtained from frequency measurements with a 
single hot-wire. The data are rather limited but sug- 
gest the following conclusions: 

(a) The scale is essentially independent of Mach 
Number (at least up to M = 0.4) and appears not to 
depend on tunnel pressure. 

(b), The intensity is greater and increases more 
rapidly with M at one atmosphere than at one-half 
atmosphere tunnel pressure. 

(c) At one-half atmosphere the intensity is low, in- 
creasing slowly with Mach Number from 0.2 per cent 
at M = 0.2 to 0.34 per cent at M = 0.7. 

The reason for (b) is not yet clear. The low values of 
intensity at one-half atmosphere pressure are gratifying 
and are somewhat surprising in view of the absence of 
damping screens in the settling region and the moderate 
contraction ratio of 8 between entrance and working 
sections. It is hoped that this extremely interesting 
turbulence investigation may be further pursued, de- 
sirable extensions including a wider variation in tunnel 
pressure, a space survey over the working-section vol- 
ume, and correlation measurements with two hot-wires. 

In connection with the turbulence studies, drag 
measurements on two spheres were made using Dry- 
den’s ‘‘pressure-sphere’’ technique.? Six- and 8.5-in. 


spheres mounted on a slender sting, stiffened by guy 


wires, were utilized. The results are plotted in Fig. 8 
for several operating pressures. Here, Ap is the differ- 
ence in pressure between orifices at the foreward stag- 
nation point and on the downstream surface of the 
sphere adjacent to the sting. The line Ap/g = 1.22 is 
indicated, since the intersection of an experimental 
curve with this line is customarily taken as defining 
the sphere critical Reynolds Number. Mach Numbers 
for several of the experimental points are noted on the 
plot. For both spheres at one atmosphere pressure the 
critical Reynolds Number lies between 350,000 and 
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360,000, which approaches the free air value of 385,000' 
but is a little lower than would be expected in view of 


the low turbulence level.* It should be noted that for 
these two cases the Mach Number is below 0.11 for all 
Reynolds Numbers up to the critical value. 

For the 6-in. sphere at one-fourth atmosphere the 
Mach Number at R = 200,000 is already nearly 0.25 
and reaches approximately 0.44 at R = 360,000. The 
pressure coefficient for this series of tests lies consist- 
ently above that for the others, and there is no sharp 
drop or any indication of a critical Reynolds Number 
as usually defined. The same phenomenon was ob- 
served and reported by Platt* (cf. Fig. 12a). It was 
there concluded that for M = 0.4 compressibility ef- 
fects radically modify the usual low-speed sphere pres- 
sure coefficient phenomena. 

In the present investigation the results with the 6- 
in. sphere at one-half atmosphere are particularly in- 
teresting in this connection. The pressure coefficient 
curve begins to drop at R = 400,000, then levels off, 
and finally has a precipitous fall at about 530,000. To 
the best of the authors’ knowledge, this is the first 
time that this phenomenon has been observed. The 
extremely high “apparent critical Reynolds Number” 
obviously*bears no relation to the usugl turbulence level 
interpretation of low-speed critical Reynolds Numbers. 
Since the Mach Number is only 0.275 at R = 385,000, 
it appears that Mach Numbers as low as this are suf- 
ficient to modify profoundly the “‘turbulence-critical 
Reynolds Number”’ relation. This is another exampl2 
of the importance of the boundary-layer-compressi- 
bility interaction which is currently under so much 
discussion. 

Summarizing the flow characteristics that have been 
discussed in this section, it appears that the CWT has 
a flow in the working section with safisfactorily uniform 
velocity distribution, the possibility of reaching Mach 
Numbers close to one, a relatively thin boundary layer, 
and a low turbulence level at high Mach Numbers. It 
should, accordingly, be particularly suitable for high- 
transonic-speed investigations. 


SPECIAL HIGH-SPEED OPERATING PROCEDURES AND 
TECHNIQUES 


The features of the tunnel drive system which permit 
air-speed control are variable-speed motors and con- 
trollable-pitch fan blade stages. The drive motors are 
a 10,000-hp. maximum output wound rotor induction 
motor and a 2,000-hp. d.c. motor with electronic speed 
control. The air speed may be controlled by a single 
toggle-switch control of power input to the motors at 
constant-pitch setting of the fan blades or by push- 
button control of the blade pitch setting at constant 


* The CWT values would appear to be in better agreement 
with the free air value of 365,000, first reported by the senior 
author,’ than with the more widely accepted value of 385,000 
given in reference 3. 


75 


2 —o6 IN. DIAM. = T 
NN | 


Bp | 4H ATMOS. 
08s | 
Lo | NUMBERS BY | 
EXPERIMENTAL POINTS, | on 
ARE MACH NUMBER | ben 
— VALUES 
| ‘ 138 | 
ee 10077 T 
| | 
ot 
2.0 3.0 40 5.0 60 


Rx 1075 


Fic. 8. Sphere pressure coefficients for two sphere sizes ai vari- 
ous tunnel pressure levels. 


shaft r.p.m. In fact, the latter method has proved so 
superior in every respect that it is now used exclusively, 
permitting exceptionally fine tolerances to be main- 
tained on speed settings and affording a considerable 
increase in operating efficiency. 

The venturi pressure difference, Ap,, which is used 
for metering the air speed, is measured by an automatic 
balance; its numerical value appears on a lamp bank in 
front of the tunnel operator in units of Ibs. per sq.ft. 
The pressure at the downstream piezometer ring, po, is 
measured on a micromanometer by a second operator. 
The ratio of Ap, to 2 has been calibrated as a function 
of the index Mach Number, M;. Direct measurement 
of this ratio with sufficient accuracy by means of a 
Machmeter is now under consideration and an appro- 
priate instrument is under construction. Meanwhile, 
a rapid and accurate graphical method for determining 
M; has served adequately for setting and maintaining 
a given Mach Number during a test run. Fig. 9 illus- 
trates the type of operating chart employed for this 
purpose. 

During vacuum operation with compressors turned 
off the general pressure level usually increases slowly 
with time because of temperature rise of the tunnel air 
and small leaks in the tunnel door and hatch seals- 
Accordingly, if constant M;, is held during a run by 
keeping both ~2 and Ap, constant, repeated and time- 
consuming operation of both compressors and speed 
control is required. If, however, the pressure is 
allowed to increase over a fairly wide range (corre- 
sponding to an approximate 5 per cent increase of 
Reynolds Number) before it is corrected by the com- 
pressors, then both Ap, and /p, can be controlled simul- 
taneously to a given ratio by the tunnel-speed control 
alone (fp, changing with Ap, as shown by the experi- 
mental curves in Fig. 9). Actual operating charts are 
much enlarged sections such as that shown within the 
square on Fig. 9. Constant Mach Number to a toler- 
ance of 0.001 is readily maintained in normal operation 
with this procedure. 
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Fic. 10. Choking condition as indicated by wall Mach Number, 
Mv, venturi pressure, Ap,, and energy ratio, E.R. 


A series of tests at various Mach Numbers can be 
run most quickly by letting the pressure vary with 
speed, using the compressors initially to establish a de- 
sired zero speed pressure and then only intermittently 
to hold the pressure within the prescribed tolerance. 
Such tests are called constant mass series, since the 
total mass of air within the tunnel shell is constant 
within tolerances dependent upon the rate of leakage 
of the tunnel seals. 

For simple modél changes and constant mass opera- 
tion 50 ten-point polars at a series of Mach Numbers 
can be run in a normal 14-hour working day. The 
running efficiency is less than this when test series are 
run at constant Reynolds Number or at constant power 
consumption. However, the loss is not great if the 
Mach Numbers are run in decreasing order when the 
tunnel is evacuated and in increasing order when it is 
pressurized, and provided that the model changes are 
long enough to permit concurrent evacuation or pres- 
surization to the minimum or maximum pressure re- 
quired. Frequently, too, the test schedule can be 
arranged so that pressurization or evacuation occurs 
during meal hours or overnight. 

The downstream piezometer ring is located a few feet 
upstream from the leading edge of the working section, 
where the cross-sectional area is a few per cent larger 
than the throat area. The area difference is sufficient 
to reduce the sensitivity of speed measurement at 
Mach Numbers above 0.9, as shown by the close spac- 
ing of the M, lines near sonic speed on Fig. 9. This has 
contributed to some of the difficulties experienced in 
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transonic testing in the normal throat. Actually, 
theory and experiment both indicate that the ring may 
be moved considerably further downstream before 
encountering model interference and that, by measuring 
the pressure on top of a small bump (approximately 
0.3-in. maximum thickness for 10-in. chord), the cali- 
bration characteristics may be considerably improved. 
Calibration with a test bump has shown that the cali- 
bration can be made four to five times more sensitive at 
M, = 1.0 by this method. 

Insistent demands for tests above 0.85 Mach Number 
in the normal throat have focused attention on the 
problems of model support, discussed in another sec- 
tion, and on the general phenomenon of tunnel choking. 
Three problems with regard to choking are of primary 
interest to the test engineer and model designer: (1) 
Although it is easy to recognize a fully choked condition 
during tunnel operation, means should be provided to 
warn of its onset; (2) it is desirable to know how near 
to choking one may approach and still obtain data that 
can be used to represent free-flight conditions; and (3) 
it would be extremely useful in model design if an 
approximate correlation could be established between 
the choking Mach Number and the geometry of the 
model and its support system. In the attempt to pro- 
vide practical solutions to these problems, wall-pressure 
data, shadowgraphs, and main-drive power data have 
been taken during tests on a great many models. 

The upper end of the experimental curve on Fig. 9 
exhibits the characteristics normally observed at 
choking. The slope of the curve, negative at lower 
speeds, becomes positive and equal to the slope of a 
constant Mach Number line, thus defining the index 
Mach Number for the fully choked condition. Further 
increases in fan power produce an increase in both 
Ap, and 2 in constant ratio, continuing until either 
the fan stalls or maximum shaft power has been 
reached. 

The change toward a positive slope, as choking is 
approached, is the result of decreasing diffuser ef- 
ficiency on the one hand and increasing air temperature 
on the other, as the tunnel seeks a higher equilibrium 
temperature for the larger input powers. When power 
input is increased with the tunnel fully choked, prac- 
tically none of the additional power is converted to 
kinetic energy and almost all appears immediately as a 
general pressure level increase. 

The experimental pressure data of Fig. 9, plotted in 
a different way, are correlated on Fig. 10 with tunnel- 
power and wall-pressure data. The particular model 
with which these data were taken had fuselage and 
nacelles mounted on an 8-ft. span unswept wing and 
was supported on a single swept strut. The axial dis- 
tribution of cross-sectional area, A, is shown for model 
and strut at the bottom of Fig. 11. It will be noted on 
the lower curve in Fig. 10 that at the lower speeds the 
energy ratio increases with Mach Number but that, as 
the speed approaches choking, the energy ratio de- 
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Fic. 11. Axial distribution of wall Mach Number, M,, at ceiling, 
and of cross-sectional area, A, of model and support strut. 


creases more and more rapidly.* The upper curves of 
this same figure show the Mach Number on the tunnel 
ceiling opposite the maximum cross-sectional area of 
model plus strut, and 1 and 2 ft. downstream from 
this position. 

Measurements show that at all points on the tunnel 
boundaries, as well as at the ceiling, in transverse planes 
just downstream of the maximum cross section, the 
air speed is supersonic when the tunnel is fully choked. 
The axial distribution of Mach Number along the ceil- 
ing is shown in Fig. 11. The speed at all points in- 
creases with M;, until sonic speed is reached opposite 
the maximum cross section. As fan power is further 
increased, speeds forward of this point remain essen- 
tially constant, while expansion to supersonic speeds 
occurs at points just downstream, followed by a com- . 
pression shock. It is now proposed that a wall-pressure 
meter, preferably a direct reading Machmeter, be in- 
stalled in the control room and connected to the wall 
orifice most suitable for use with a given model installa- 
tion. This instrument would then provide the neces- 
sary warning of the onset of choking with a sensitivity 
not inherent in the other available methods (Ap, and 
energy ratio). The further improvement in tunnel 
operation to be gained from the use of this technique is 
discussed in a later section. 

The data of Figs. 10 and 11 indicate that the tunnel 
speed may be brought to within 0.01 of the choking 
Mach Number before sonic speed is reached at the 
wall. The same result has been found for all models for 


* For this large model a maximum energy ratio of 7.0 was ob- 
tained as compared with 7.5 normally measured for a clear throat. 
The energy ratio values are based upon shaft power. 
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Fic. 12. Choking Mach Number as function of blocked area. 


which pressure data are available; this includes models 
giving choking Mach Numbers up to 0.95. Shadow- 
graph pictures have shown the rapid growth of the 
supersonic region about a model near choking speeds. 
Shock waves at the model extend rapidly toward the 
walls beginning with Mach Numbers approximately 
0.01 less than choking M;. 

The correlation between choking M; and maximum 
blocked area given on Fig. 12 includes data for a great 
many different model configurations and supports. 
The experimental data are divided into five major 
groupings and are compared with the theoretical curve 
of choking Mach Number computed from one-dimen- 
sional theory. In the clear tunnel, choking occurs at 
-the downstream end of the working section, limiting 
the value of 17; to 0.98 in two of the working-section 


Fic. 18. Lockheed bump mounted in the CWT with velocity 
survey plate in model position. 


carts and to 0.96 in the third cart, as indicated on the 
figure. The blocked area is defined as the maximum 
cross-sectional area taken from a distribution such as 
that shown at the bottom of Fig. 11. The distribution 
from the model is added to that from the swept strut 
or struts in exact accordance with the position of each 
along the axis of the tunnel, and the maximum value 
of the sum is taken as the blocked area. The experi- 
mental points are in fair agreement with the theoretical 
curve. The points deviating furthest from the curve 
are those for a swept strut tested with small tip fairings 
and for medium-sized bodies of revolution tested on 
this same strut. With these configurations it was pos- 
sible to exceed the theoretical choking limit by about 
0.03; with all other configurations, the tunnel choked 
within 0.01 of the theoretical M;. 


THE LOCKHEED TRANSONIC BUMP 


The Lockheed “bump” has proved of notable value 
for transonic testing in the wind tunnel and deserves 
special consideration. Complete details of its charac- 
teristics and operation have been presented elsewhere.° 
Briefly, the bump is a cambered profile mounted on the 
floor of the test section and extending from wall to 
wall (Figs. 13, 14). Between the floor and the under- 
surface of the bump is a boundary-layer removal duct, 
installed to by-pass the comparatively thick boundary 
layer developed along the tunnel floor ahead of the 
bump; the boundary layer on the upper surface of the 
bump at model position is approximately 1/2 in. thick. 

Small semispan models mounted with span normal 
to the bump surface (reflection plane mounting) are 
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tested in the high-velucity field developed locally in the 
vicinity of the bump maximum thickness. The tech- 
nique is closely related to the N.A.C.A. wing-flow 
method used in flight, from which published data are 
already available. Both methods take advantage of 
the local regions of transonic or supersonic velocity, 
supported by the flow past curved surfaces, in which 
small models may be tested with continuously variable 
speed throughout the transonic range. In comparison 
with the wing-flow technique the bump method has all 
of the advantages associated with more conventional 
wind-tunnel methods: flexibility in and accurate con- 
trol of operating conditions, relatively easy instrumen- 
tation, economy and rapidity in obtaining data, etc. 


The speed range available for tests on the bump is 
limited at the high end by tunnel choking but extends 
to a maximum Mach Number of about 1.15, for which 
the Mach Number just ahead of the bump is a little 
over 0.6. Higher performance might be obtained from 
a bump of smaller chord and greater camber; however, 
such a configuration would almost certainly be limited 
by shocks at the model position if not by choking. The 
steeper velocity gradients associated with greater 
camber would also reduce the size of the models which 
could usefully be tested. 


The average wing span of the models that have been 
tested with the Lockheed bump is approximately 6 in. 
Strong shocks occur on the bump, but these are down- 
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Fic. 14. Diagram of Lockheed bump, and Mach Number dis 


tribution over region occupied by model. 


stream of the model position and have produced no 
noticeable additional unsteadiness in the flow at the 
model. Actually, the entire flow through the working 
section is subject to high-frequency velocity fluctua- 
tions of approximately '/. per cent amplitude because 
of the large diffuser losses induced by the bump and 
the separation near the bump trailing edge. 
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Fic. 15. Drag coefficient for three similar swept wing models in the transonic region, using different testing techniques. 
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The only limitation on the low end of the speed range 
is imposed by the reduction in magnitude of the forces 
acting on the models, rendering their measurement more 
difficult. In a variable-density wind tunnel this may 
be avoided by operating at a constant dynamic pressure 
for all Mach Numbers, provided that the resultant 
variations in Reynolds Number are acceptable. The 
Reynolds Number variation with Mach Number for 
constant dynamic pressure is opposite to, and consider- 
ably greater than, that for the constant mass operation 
discussed earlier. 


The velocity field in the vicinity of a model at a 
nominal Mach Number of 1.04 is shown in Fig. 14, 
where a typical wing model is also shown mounted on 
the bump. The decrease of velocity from root to tip of 
the sample model is a little over 2 per cent; over a 
10-in. semispan the decrease would be 5 per cent. 
This lateral gradient is essentially the same for all 
speeds in the transonic range. In the direction of flow 
there are small deviations in velocity; the Mach Num- 
ber has 1 per cent less than its centerline value at dis- 
tances 3 in. forward and 5 in. aft of the model center- 
line. The longitudinal gradients become somewhat 
larger at higher supersonic speeds but become even less 
at subsonic speeds. Test results are referred to an 
average speed or nominal Mach Number taken in the 
region occupied by the model (Fig. 14). Approximate 
allowance for the '/2-in. thick boundary layer can be 
made by building 0.2-in. thick inboard extensions on 
the models as indicated on the figure. Tests need not 
be limited to wing models, since semifuselages and semi- 
span horizontal tails can be included as well, thus per- 
mitting studies of wing-fuselage interference and tail 
trimming effects. The Reynolds Numbers for bump 
tests, like those in the wing flow method, are rather 
small (from 500,000 to 1,000,000 at M@ = 1.0). How- 
ever, this should not detract from the great usefulness 
of the bump as a means of exploring the transonic zone. 


Measurements of forces and moments on the models 
are made with a small six-component strain gage bal- 
ance with self-contained pitch-changing mechanism; 
the balance rotates with the model and thus measures 
the forces relative to body axes. Typical drag results 
for a swept-wing model are shown on Fig. 15, giving the 
drag at zero lift as a function of Mach Number through 
the transonic range. These measurements were taken 
at constant dynamic pressure on a 5-in. semispan model 
and indicate the quality of results yielded by the balance 
system. Lift and pitching moment results are equally 
good; the data in general are repeatable to better than 
1 per cent. The models can be constructed quickly and 
cheaply and can readily be installed or removed. This 
feature, together with the high performance of the bal- 


. ance system and tunnel speed control, has permitted 
_the running of 65 seven-point polars at various Mach 


Numbers and constant dynamic pressure in a normal 
14-hour working day. 


The other two curves shown on Fig. 15 give the re- 
sults of tests on larger models of the same wing run in 
the normal subsonic working section. The lower curve 
gives the drag measured on a 2.5-ft. semispan model 
mounted on the floor turntable of the working section, 
while the central curve gives the drag measured on a 
complete 5-ft. span wing mounted on a sting extending 
forward from swept struts. Blocking corrections have 
been applied to these data according to the method out- 
lined in the final section of the paper. The sting- 
mounted model could not be carried quite up to the 
critical Mach Number, but the bump and reflection- 
plane tests show drag rise at very nearly the same Mach 
Number. The difference in drag level between the 
bump and the sting tests is largely explained by the dif- 
ference in test Reynolds Numbers. The tares for the 
sting and struts were accurately determined, whereas 
those for the reflection-plane mount were not. How- 
ever, it is known that the reflection-plane tares are 
nearly independent of Mach Number. Hence, the drag 
curve shape for the reflection-plane test should be very 
nearly correct, any tare errors having only the effect of 
shifting the curve vertically. The similarity in the 
three sets of data near and below the critical Mach 
Number, where the drag begins to rise sharply, suggests 
that considerable confidence may be placed in the bump- 
test data at the higher Mach Numbers, which the other 
two testing techniques cannot reach. 


METHODS OF MODEL SUPPORT 


One of the major problems confronting the wind- 
tunnel operator is that of devising a method of support- 
ing the model in the wind stream under severe air loads 
while at the same time keeping the supporting structure 
small enough to permit accurate determination of the 
model characteristics. Even for low-speed operation 
the support problem is a difficult one, requiring lengthy 
and frequent tare tests. For high speeds where com- 
pressibility effects are pronounced, the model supports 
must be designed to eliminate as far as possible local 
high-speed flow regions with their accompanying shock 
waves. Both the model and its supports must be small 
enough to allow the attainment of the desired M; before 
tunnel choking occurs. Conversely, it is desirable to 
keep the model size as large as possible so as to avoid 
operating near or below the critical Reynolds-Number 
range, and to generate air loads of a magnitude that 
can accurately be read with the available balance sys- 
tem. Also, the model-supporting system must have 
sufficient strength and rigidity to carry these air loads 
without failure or excessive deflection. Clearly, the 
aerodynamic and structural requirements are in direct 
conflict, as is usually the case in aeronautics, and the 
wind-tunnel operator must resort to compromises and 
“trick’”’ solutions. 


The most important conditions that a model and 
support system must satisfy are: (1) sufficient strength 
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(2) small enough blocking area so that 


and rigidity; 
choking will not occur in the test Mach Number range 
desired; (3) tares negligibly small, accurately measur- 


able, or essentially constant; (4) no critical interfer- 
ence effects; (5) adaptability to the available balance 
system if force and moment data are required. 


In the CWT, condition (1) is based on the require- 
ment that neither the model nor the support system 
should fail structurally. Some German tunnels have 
“model catching screens” installed downstream of the 
tunnel throat to protect the tunnel structure and fans 
in case of failure. This permits a considerable increase 
in the load limits on the support system, since it is ex- 
pected that some fraction of the models will be lost be- 
cause of support-system failure. Under such condi- 
tions the design strength is primarily determined by 
economic considerations. The CWT does not subscribe 
to this philosophy and, accordingly, uses a safety factor 
of five based on the ultimate strength of the material. 


Condition (2) must be satisfied since the experi- 
mental data are meaningless if the tunnel is choked. 
Conditions (3) and (4) refer to the direct air loads on 
the support system and to the interference of the latter 
on the model air loads. Both of these effects must be 
determined and subtracted from the total measured 
loads. Great care must be used to avoid configurations 
where the support system appreciably influences the 
occurrence of such critical model phenomena as bound- 
ary-layer transition, shock waves, and flow separation. 
Condition (5) is simply the requirement that the model- 
support system must be mounted on one of the avail- 
able balance systems. This condition frequently rules 
out the use of many types of model support. 

The five conditions given above have been set up, 
based on the assumption that the objective of the test 
is to obtain accurate force and moment data; this is 
by far the most common type of test. For other types 
of tests, such as those with pressure and duct models, 
condition (3) may not be important. Finally, consider- 
ation should be given to the cost of a proposed support 
relative to the value of the results, to the use of free 
operating parameters such as angle of attack and yaw, 
and to applicability to more than one model. 


The most important types of model-support system 
are: (1) cantilever strut mountings, (2) wing tip and 
reflection plane mountings, (3) fixed-end beam and 
center-plane mountings, (4) sting mountings, and (5) 
wire and tension strap mountings. 

The, cantilever strut mountings are probably the 
most commonly used support systems today, at least 
for three-dimensional tests. The models are attached 
to the free ends of the struts, and the strut bases are 
attached to the balance system. Either one, two, or 
three struts may be used, all being attached to the 
lower side of the model. For low-speed tests and for 


tests involving large cross-wind forces and yawing 
moments, the CWT uses wind-shielded vertical struts 


(strut axis perpentlicular to the flow direction) with 
moderate fineness ratios. For high-speed tests un- 
shielded swept struts are required. The CWT has 
three sets of swept cantilever struts: (a) a three-strut 
arrangement with image system, (b) a single strut with 
image system, and (c) a single strut with an internal 
pitch changing link. All struts have sweep angles be- 
tween 35° and 40° and thickness ratios between 10 
and 12 percent. Apparently reliable results have been 
obtained with these swept struts up to If; = 0.90. 
It is doubtful if they can be used for M, > 0.93 unless 
their blocking effects are nullified by modifications to 
the cross-sectional area of the working section. 

Experience at the CWT has shown that vibration 
of both model and strut is usually the limiting struc- 
tural condition, particularly for the single strut mounts. 
Great care is used to avoid the occurrence of large os- 
cillations that might result in elastic instability,* and 
all force and moment balances are placarded before 
the start of a test Examples of the strut tares and 
interference effects are given in the next section. 

A large portion of the 30 tests run at the CWT have 
made use of the wing-tip and reflection-plane support 
systems. For this type of support the model projects 
through the tunnel wall and attaches directly to the 
balance system, thereby eliminating the strut tare, 
interference, and strength problems. Strength of the 
model is not a serious problem for reflection-plane 
models (where the tunnel wall forms the center plane 
of the model) but is an.important restriction for wing- 
tip mounts. The chief experimental difficulties arise 
from the effects of the tunnel-wall boundary layer,® 
of the air loads on the turntable fairing covering the 
gap between model and tunnel wall, and of the seal 
between the model and the tunnel wall or fairing. No 
attempt has yet been made at the CWT to measure 
these effects directly, even though low-speed tests in 
the GALCIT Ten Foot Tunnel have shown the drag 
tare to be of considerable magnitude, especially for high 
lift conditions. For reflection-plane setups with a 
wing-fuselage-tail configuration and high Mach Num- 
ber, it is expected that the boundary layer will have a 
large effect because the model must be small to avoid 
tunnel choking. In some recent but unfinished tests 
at the CWT, boundary-layer removal devices are being 
utilized and indicate promising results. 

Several tests have been made with the model mounted 
on a swept-wing strut that completely spanned the 
jet and attached to the balance system outside of the 
tunnel walls. The model was mounted at the center 
of this ‘‘fixed-end beam.’’ This system has greater 
strength and rigidity than does a cantilever strut ar- 
rangement but in some cases may have a greater inter- 
ference effect. For center-plane support systems it is 
expected that the interference and wall boundary-layer 

* The lateral buckling equation for these doubly tapered struts 


has been solved mathematically by CWT personnel. The results 
will be published in the near future. 
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Fic. 16. GALCIT image system tare procedure. Step 1. 
Make polar Runs 1, 2, 3, and 4, covering full angle of attack 
range at constant free-stream dynamic pressure. Step 2. Cor- 
rect all data of Runs 2 and 3 for average increase in dynamic pres- 
sure due to presence of image system as determined by difference 
in lift curve slope between Runs 1 and 2, and 3 and 4. Step 3. 
Subtract results of Run 1 from averaged results of Runs 2 and 3. 
Step 4. Values of Step 3 are the tares for model normal and are 
to be subtracted from the model normal data of Run 4. 


effects will be small but that the tare problems will be 
the same as for a reflection-plane system. The center- 
plane support has an advantage over the reflection- 


plane support only insofar as the wall boundary-layer | 


effects are concerned. 


Sting mountings cause little interference and block- 
ing (at the model) and have nearly constant tares. 
The blocking of the sting support can be largely nullified 
by proper modification of the tunnel cross-sectional 
area. These important advantages are offset by an in- 
herent lack of strength and rigidity, which can, how- 
ever, be considerably improved by the use of stabilizing 
wires. The CWT is currently constructing a sting 
mounted on a swept strut which will utilize a three- 
component, strain gage balance located inside of the 
model. 


The wire or tension strap support system would ap- 
pear to have many advantages and might well satisfy 
the first four conditions better than any of the other 
systems mentioned. Its most serious difficulty is that 
it requires a ring or equivalent structure mounted on the 
balances and completely surrounding the working sec- 
tion. This ring structure must also extend a consider- 
able distance up and downstream from the model to 
permit the attachment of support wires with large 
sweep. The CWT ring frame support does surround 
the working section but does not extend far enough up 
and downstream. It is hoped that it may be modi- 
fied so as to make possible tension strap support 
tests. 


The preceding discussion has been concerned pri- 
marily with support systems used or proposed for tests 
requiring force and moment measurements. Special 
systems devised when necessary for unusual test con- 
ditions are nearly always modifications of the types 
mentioned above. It has become evident that there 
has been no support system developed up to this time 
which is usable for more than a few model types and 
test conditions. For technically and economically 
satisfactory performance the transonic wind tunnel 
must build up a complete stockpile of different support 
systems with their associated balances. This is a long- 
range and expensive project, which can only be carried 
out over an extended period of time. 


IMAGE SYSTEM TARE TESTS WITH CANTILEVER STRUT 
SUPPORTS 


For cantilever strut supports the most reliable tare 
procedure seems to be the image system method. Un- 
fortunately, this method involves making three to five 
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Fic. 17. Swept strut image system tare procedure. Step 1. 
Make polar Runs 1 through 6, covering full angle of attack 
range at constant M;. Step 2. Add results of Runs 2, 3, 5, and 6 
together at constant angle of attack values. Step 3. From half 
of the summation of Step 2 subtract the results of Run 1. 
Step 4. Values of Step 3 are the tares for model normal and 
are to be subtracted from the model normal data of Run 4. 
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extra runs for every major model change and therefore 
is expensive. A thorough image system tare program 
has been carried out at the CWT for only one model so 
far. Examples of the results of these tests are described 
below. 

The model was mounted on two different cantilever 
strut systems: (a) three-member vertical struts with 
windshields, for low Mach Number tests; and (b) 
three-member swept struts without windshields. For 
the vertical struts the GALCIT tare procedure was 
used. The method is briefly outlined in Fig. 16. The 
image system duplicates the main support struts and 
windshields in size, shape, and relative location. The 
image bayonets (removable tips of the struts) are fast- 
ened to the main bayonets and always remain clear of 
the image windshields, so that all air loads on the image 
bayonets are transmitted to the balance system. The 
basic principle (and assumption) of the tare procedure 
is that the addition of the image system with inverted 
model changes the balance readings by the amount of 
the tares of the main struts with model normal. To get 
the free-stream inclination and curvature correctiops, 
the runs with image system in, model normal and in- 
verted, are compared. All data are compared at the 
same values of the lift coefficient. Drag and pitching 
moment tares are shown on Figs. 18 and 19. The lift 
tares are small (C,;, < 0.06). 


For the swept struts a modification of the GALCIT 
tare procedure was used. The method is briefly out- 
lined in Fig. 17. In this figure the tips of the image 
struts are depicted as if they did not extend to the 
model, to show that they were not in physical contact 
with the model. Actually, the tips extended inside 
the model at the trunnion wells. The procedure cor- 
rects for free-stream inclination and curvature effects 
by a method simila: to that used for the vertical strut 
tares. With these particular swept struts and model 
there is no measurable effect on the average dynamic 
pressure, as indicated by the slope of the curve of lift 


versus angle of attack. Thus, all data can be com-. 


pared at the same values of angle of attack. For com- 
parison, forces and moments were also measured on the 
swept struts alone without any model in the tunnel. 
Drag and pitching moment tares are shown on Figs. 
18 and 19. 


The drag tare curves of Fig. 18 indicate several inter- 
esting characteristics. First, the inclination effect of 
the vertical struts and windshields, as shown by the 
negative slope of the tare curve, is not present in the 
swept-strut tares. Second, the vertical strut tare at 
C, = 0 is higher at low Mach Numbers than the corre- 
sponding tare for the swept struts. Third, for low Mach 
Numbers the drag of the strut alone is approximately 
equal to the measured tare, but for high Mach Numbers 
and angles of attack the drag interference between the 
swept struts and the model makes up a large percentage 
of the total tare. 
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Fic. 20. Mach Number variation of swept strut tares. 


The pitching moment tare curves of Fig. 19 show 
tares of opposite signs for the two types of strut. This 
probably indicates an inclination effect with the vertical 
strut system, for which the inclination of the wind 
stream caused by the large windshields changes the 
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Fic. 19. Pitching moment tares for vertical and swept struts. es 
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Fic. 21. Comparison of flight-test and wind-tunnel results for 
a bomber aircraft. 


model tail load. The “swept struts alone” tares are 
appreciably different from the image system tares. 

In Fig. 20 are shown, for two angles of attack, the 
variations with Mach Number of the drag and pitching 
moment tarés as determined by the “struts alone’ and 
image methods. It is clear that a determination of 
tares by the simple procedure of force and moment 
measurements on the support system without model 
present may lead to considerable errors, particularly for 
drag at high Mach Numbers. 


COMPARISON OF WIND-TUNNEL RESULTS WITH FLIGHT 
TESTS 


The ultimate goal of most wind-tunnel tests is to 
enable the airplane designer to predict the performance 
of his airplane in actual flight. To get a truly reliable 
comparison is most difficult because the wind-tunnel 
model and the airplane seldom have exact aerody- 
namic similarity. The surface condition of the airplane 
is almost always rougher than that of the model. The 
airplane has antennae, pitot-static tube mounts, drain 
pipes, air vents, flexible skin covering, air leakage, and 
a great many small protuberances, all of which add up 
to an appreciable portion of the airplane drag. Power- 
plant performance is difficult to measure precisely. 
Characteristics of control surfaces are often sensitive to 
minor changes in contours and gap shape and size. In 
practice, the majority of wind-tunnel tests are set up to 
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determine the comparative effects of various modifica- 
tions which the airplane designer then correlates with 
his theoretical calculations and such other experimental 
data as are available to build up his performance 
predictions. A few serious attempts have been made 
to obtain a direct correlation of flight and wind-tunnel 
data.’ Some recent flight-test results have been 
gathered by the CWT for comparison with other tests 
in both the CWT and GALCIT Ten Foot wind tunnels. 
These data are shown in Figs. 21-23. In all cases the 
origins have been shifted to obscure the true level of the 
data, but differences between the points as shown and 
all of the Mach and Reynolds Number values are cor- 
rect. 

For the wind-tunnel tests of Fig. 21 only strut-alone 
tares were measured, since there was no image strut 
available. Comparison with low-speed tests of the 
same model in another tunnel, wherein complete image 
system tares were determined, shows that 90 per cent 
of the differences between tunnel and flight lift and 
trim data are caused by the incomplete tare determina- 
tion in the CWT tests. For the aileron effectiveness 
plot, the rolling resistance data of Pearson® were used to 
calculate the roll parameter, pb/2V, from the tunnel 
data. 

For the wind-tunnel tests shown in Figs. 22 and 23, 
complete image system tares were measured. The 
flight-test drag data were corrected for the drag due to 
engine cooling, ventilating system, antennae, and mis- 
cellaneous protuberances. The agreement between 
wind-tunnel and flight data is unusually good. 


WIND-TUNNEL WALL EFFECTS IN INCOMPRESSIBLE 
FLow 


The boundary effects produced by wind-tunnel walls 
are primarily the increase in axial velocity at the model 
due to constriction of the flow (blocking) and the 
boundary induced up-wash due to lift. Solutions for 
these effects in incompressible flow have been obtained 


O CWT, M=0.300,R=1.3 x 108 | 
CWT, M=0.203,R=3.8Xx 10 


— — FLIGHT TEST, M = 0.270, R=26.0 x10® 
i 
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Fic. 22. Comparison of flight-test and wind-tunnel results for a 
transport aircraft. 
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by many investigators, and some of these are applicable 
directly to the geometry of the CWT working section. 
The latter consists of a basic rectangle, whose sides are 
in the ratio 1/2, converted to an octagonal section by 
the addition of rather large corner fillets. Calculations 
are most easily carried out first for the basic rectangle, 
permitting solution by a doubly infinite system of 
image vortices for the up-wash corrections or by a sys- 
tem of image doublets for the blockage corrections. 
Small modifications to these results are then made to 
account for the fillets. 

Methods outlined by Batchelor® have been applied to 
the problem of up-wash corrections for straight rec- 
tangular wings, leading to numerical factors for use in 
the following formulas for the up-wash corrections to 
drag, angle of attack, and pitching moment. 


ACp = 6(S/BH)C,? 
Aq = [1+ (1) 
AC, = 1(S/BH)(dC,,/ds)C, 


Here, B and H/ are the tunnel breadth and height, re- 
spectively; Sis the reference area for the lift coefficient; 
and cis the MAC of the wing. The other notations 
have their usual meanings. The wind tunnel for which 
Batchelor has carried out calculations for the effect of 
corner fillets is almost identical in geometry of cross 
section to the CWT. It is possible to use his results 
directly, with slight modifications to account for the 
difference between corner fillets. The dimensionless 
factors 6, k, and r are evaluated for the CWT geometry 
for various sized models. The numerical factors for 
centrally located models are also applicable to re- 
flection-plane models mounted on the floor or ceiling 
turntables of the CWT wing-tip cart. 

Up-wash corrections for sweptback wings have been 
investigated” for a circular tunnel. Corrections for 
yawed wings in a 7 by 10 rectangular tunnel are given 
by Swanson.!! Similar calculations for arbitrary span- 
wise and chordwise loadings of sweptback wings are 
now being carried out for the CWT. 

Lock” has derived the solid blockage corrections for 
bodies of revolution and for two-dimensional wings, in 
rectangular tunnels. Numerically identical results are 
derived by Thom,'* who introduces model volume into 
the formulas and who also carried out calculations for 
finite span rectangular wings. The results for bodies of 
revolution and two- and three-dimensional wings can be 
averaged without much loss in accuracy,'* so that a 
single formula may be used for all three, as well as for 
complete airplane models. The formula expressed in 
terms of model volume, and with the average numerical 
factor that has been found to fit most caseg encountered 
at the CWT, is 


= (Au/U)sotia = 0.89 [Volume/(BH)"] (2) 


Calculations for the blockage of swept wings have so 
far not been published. These are being carried out 
at the CWT at the present time, as are also the calcu- 
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Fic. 23. Variation with Reynolds Number of wind-tunnel model 
parasite drag and maximum lift, and comparison with flight test. 


lations for the velocity at the tunnel wall opposite the 
model. Wall velocity calculations have been given by 
Thom" for straight wings and have been extended by 
other authors. Correlation between the velocity at the 
wall and the central blockage is needed for the high 
Mach Number wall correction procedures discussed in 
the final section of the paper. 
The formula for the blockage of the model wake*’ 1° 
is 
= (Au/U) ware = '/4(SCpp/A) (3) 


where S is the reference area for the drag coefficient, 
A is the cross-sectional area of the octagonal tunnel 
throat, and Cp, is the measured parasite drag coefficient. 
Another correction for the wake effect is the buoyancy 
correction to drag due to the pressure gradient induced 
at the model by the wake.” Heretofore considered 
negligible, its effects are now being computed at the 
CWT because of its relative importance at extremely 
high subsonic Mach Numbers. 

If not evaluated by tare measurements as outlined 
in a previous section, the blockage and interference ef- 
fects of the model-support system must be computed 
from similar considerations. 


EXTENSION OF WALL EFFECTS TO MODERATE SUBSONIC 
Macu NuMBERS 


The procedures outlined above for the determination 
of wind-tunnel boundary interference effects are valid 
as long as all velocities are small in comparison with the 
speed of sound. As tunnel speeds increase, compressi- 
bility effects enter which modify, at first slightly and 
later profoundly, these tunnel-wall corrections. The 
only theoretical treatment that has as yet been at all 
fully developed in this connection is that involving 
linearization of the differential equations by the method 
of small perturbations. The so-called Prandtl-Glauert 
rule was developed in this way and has been widely 
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used. A complete treatment, with particular applica- 
tion to wind-tunnel interference, was given by Gold- 
stein and Young” and has been used as a basis for many 
subsequent studies.'* Unfortunately, there has been a 
good deal of confusion in connection with the applica- 
tion of the Prandtl-Glauert rule to certain special cases, 
especially to the flow about slender bodies of revolu- 
tion.” In fact, the compressibility effect on certain of 
the wind-tunnel corrections for such bodies is given in- 
correctly by Thom" and Goldstein and Young.'® 

In view of this confusion the authors have found 
particularly clarifying and useful a formulation given 
by Goethert which he calls the “extended Prandtl 
rule.”’* As applied to wind-tunnel configurations this 
rule may be stated as follows: 

“Given a model in a wind tunnel with compressible 
fluid flow. Let the x axis be in the wind direction with 
the y and z axes perpendicular. The velocity far ahead 
of the model is taken as U, the Mach Number as M, and 
B is defined as V1 — M?. Construct a comparison 
system with incompressible flow by decreasing all y 
and zg dimensions of model and tunnel by the multi- 
plicative factor 8 while leaving all x dimensions un- 
altered. This reduces tunnel diameter, model span, 
model thickness, and angles of attack and yaw by the 
factor 8, while leaving model chord or body length un- 
altered. Calculate normal tunnel-wall corrections for 
the comparison system using incompressible fluid 
theory. These corrections have the form Au/U, 
Av/U, and Aw/U corresponding to axial, side-wash, 
and down-wash velocities. Obtain the compressible 
flow corrections by multiplying the incompressible 
values at corresponding points by the following factors: 
Au/U by 1/8, Av/U and Aw/U by 1/8.” 

This rule makes it possible to determine easily and 
without chance of confusion the compressible flow 
tunnel-wall corrections for any case where incompres- 
sible flow calculations have been carried out. Applying 
the rule to the correction formulas collected above leads 
to the ‘‘compressible”’ correction formulas used at the 
CWT: 


AC), and Aa of Eqs. (1) are unaltered by com- 
pressibility* 

AC,, of Eqs. (1) is altered only in that the coef- (4) 
ficient 7 is determined for a tail length longer 
than the actual one by the factor 1/8 


(exe = (€s)ine./B* (5) 


where the subscript c indicates the actual correction 
used for compressible flow and inc. denotes the correc- 
tion calculated from the true model and tunnel geome- 
try for the case of incompressible flow, using Eq. (2) 


= (1/B*)(*/4)(S/A) (Cop). (6) 


* In the equation for Aa, the term with the factor k is a small 
correction for the induced flow curvature over the wing. In 
compressible flow, for small c/H, the factor k is increased by 1/8. 
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where (Cp,), is the measured drag coefficient in the 
compressible flow. 


TUNNEL BLOCKAGE EFFECTS FOR HIGH SUBSONIC 
Macu NuMBERS 


Up to tunnel Mach Numbers of the order of 0.8, the 
above blockage correction formulas [Eqs. (5) and (6)] 
are satisfactory for slender, well-streamlined models 
and support systems whose dimensions are moderate 
compared to the tunnel diameter. Unfortunately, as 
M approaches 1 (6 — 0), the linearized theory is no 
longer valid, as is made obvious by the fact that the 
“small’”’ correction terms approach infinity. There is 
as yet no transonic theory far enough developed to be 
of any assistance in this connection, so that the wind- 


tunnel operator is forced for the present to combine | 


empirical methods with such extrapolations of the linear 
theory as appear not too unreasonable. Such a com- 
bination that seems promising is afforded by the use 
of measured “wall Mach Number,” //,, as basic in- 
dependent or operating parameter, rather than index 
Mach Number, M;,, as determined from piezometer 
pressures far upstream from the model. 

A practical calculation difficulty appears if the at- 
tempt is made to compute free-stream Mach Number 
at the model by applying the linearized theory blocking 
corrections to the measured M, under conditions ap- 
proaching that of tunnel choking. The 8 which appears 
in the denominator of the correction formulas (5) and 
(6) should, even as a first approximation, correspond 
to some sort of average Mach Number over the trans- 
verse plane containing the maximum cross-sectional 
area of the model. An iteration procedure must ac- 
cordingly be followed. 

It is found in carrying out such computations that, as 
M, increases, a point is reached where an infinitesimal 
increment in M;, produces a sudden jump in the cal- 
culated free-stream Mach Number at the model, from a 
value well below M = 1 to one that is definitely super- 
sonic. This computational difficulty is related to an 
experimental observation made by Dr. Liepmann dur- 
ing recent fundamental studies of transonic-flow 
phenomena at the California Institute Guggenheim 
Laboratory. He found that it was practically impos- 
sible to duplicate conditions over the surface of an air- 
foil by setting at a definite value the Mach Number in 
the undisturbed flow far ahead. However, if the Mach 
Number were set at a fixed point in the local high veloc- 
ity region near the airfoil, the entire flow pattern could 
accurately be reproduced. This means that in the 
transonic regign the flow conditions are enormously 
sensitive to small changes in the undisturbed stream 
Mach Number and, from a computational point of 
view, that calculations of a higher Mach Number from 
an observed lower one are highly untrustworthy. On 
the other hand, the calculation of a lower Mach Num- 
ber from an observed higher one should be relatively 
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safe. Theory and experiment both indicate that the 
corrected ‘‘free-stream’’ Mach Number, M, at a model 
in a closed wind tunnel is higher than the Mach Num- 
ber, /,, in the undisturbed flow far upstream but lower 
than the Mach Number, M,, at the wall in the trans- 
verse plane containing the model. 

As has already been mentioned, Thom has discussed 
the use of wall pressures, and other British investigators 
have subsequently studied the matter further. Earlier, 
Goethert” showed that the extended Prandtl rule could 
legitimately be used to calculate blocking corrections 
from observed wall pressures, even though the linear- 
ized theory assumptions did not hold in the neighbor- 
hood of the model. The essence of Goethert’s argu- 
ment is that the ratio, 1/m, of the velocity increment at 
the model due to blocking, to the disturbance velocity 
produced by the model at the wall is, to a first approxi- 
mation, independent of Mach Number. Furthermore, 
even in the extreme cases where the variation of this 
ratio with Mach Number is important, the extended 
Prandtl rule gives a method of obtaining a good esti- 
mate of the effect. 

The above considerations, combined with extensive 
experimental studies made during the course of many 
routine tests and special calibration investigations, have 
led to the following procedure for obtaining blocking 
corrections in the CWT at high subsonic speeds (the 
procedure is outlined for models at zero lift; the exten- 
sion to cases with lift involves only the averaging of 
pressures from orifices at two opposite walls) : 

(1) Using pressure orifices on the wall opposite the 
maximum cross-sectional area of the model and near 
the model plane of symmetry, determine empirically the 
pressure difference, Ap,, produced by the model. Plot 
Ap,,/ Ap, vs. Mi. 

(2) Compute the wall pressure increment due to 
wake blockage (constant across the tunnel cross sec- 
tion) from the measured parasite drag and the wake 
blockage formula [Eqs. (3) or (6)]:* 


Ap(wake) /q; 2(Ew)e 
Plot 


Ap(wake) Ap(wake) ( ) 
Ap, Ap, 

(Here g; is the index dynamic pressure, determined as 
a function of Ap, and M;, by the clear tunnel calibra- 
tion. 8 is taken as V1 — My?.) 

(3) Determine the portion of Ap, arising from solid 
blocking: 

Ap,..(solid)/ Ap, = (Apw/ Ape) — [Ap(wake)/Ap,] 

(4) Determine the solid blockage correction [corre- 
sponding to eqs. (2) and (5)] from the corrected wall 
pressure: 


* This is the only step in the procedure where linear theory 
calculations are fully used. The factor here computed is nor- 
mally small, and divergence in the calculations is avoided by the 
definition of 


2Ap, 


where m is the ratio of disturbance velocity at the wall 
to the solid blockage velocity increment at the model, 
calculated on the basis of incompressible fluid. Numer- 
ical values for a number of configurations are given in 
references 13 and 19. Other configurations including 
swept wings and struts are now being calculated at the 
CWT. A good average value for m is 3. 
(5) Determine the total blockage correction 


€ = + 
and plot e, vs. M;. 
_ (6) Finally, compute the corrected “‘free-stream” 
quantities 
M = M,{1 + «(1 + 0.2M,’)] 
gq = + «(2 — M,’)] 


The factors in parentheses take account of the adia- 
batic expansion resulting from the blocking velocity 
increment. 

(7) The “‘free-stream” coefficients are obtained by 
dividing measured forces and moments by g (and the 
appropriate area or volume) and should be plotted 
against the ‘‘free-stream’’ Mach Number, M. 


—— FROM WALL PRESSURES 
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Fic. 24. Wind-tunnel blocking correction factors for a particular 
model by the wall-pressure and linearized theory methods. 


Fig. 24 shows the results of such a data reduction 
procedure in the case of the semispan, reflection-plane 
wing model of Fig. 15. The corrections computed on 
the basis of the linearized theory are also presented. 
For this configuration, with a small model and no sup- 
port blockage, the corrections remain moderate up to 
the actual choking point. Note that the corrected free- 
stream Mach Numbers as indicated by the linearized 
theory are less than those given by the wall-pressure 
method. Thom" has already suggested that the linear 
theory probably gives too small free-stream Mach 
Numbers, and the present results support this conjec- 
ture. For larger models and strut support systems both 
the correction factors and the differences between the 
two procedures become rapidly greater. 


the 
(6)] 
dels 
rate 
as 
no 
the 
is 
o be Bie: 
ind. 
bine 
near 
in- 
idex 
eter 
 at- 
aber 
king 
ap- 
ears 
and 
yond 
ans- 
onal 
t, as 
imal 
cal- 
ym a 
dur- 
flow 
leim 
pos- 
air- 
in 
[ach 
; 
ould 
the 
usly 
eam 
t of 
rom — 
On 
um- = 
vely 


For strut-mounted model tests wall pressures must be 
observed both with support struts alone and with 
struts plus model. The data correction procedures are 
somewhat more elaborate but are identical in principle 
with those outlined above. The data of Fig. 15 for the 
wing so mounted have been corrected in this way. 


CONCLUSION 


The present paper makes it apparent that much 
further research, both experimental and theoretical, is 
required before transonic phenomena, or even the inter- 
pretation of transonic wind-tunnel tests, can be properly 


understood or data used with full confidence. It is 


hoped that the techniques and procedures here pre- 
sented may be found useful in this further development. 
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N.A.C.A. Cascade 
Design of High-Performance 
Axial-Flow Compressors 


SEYMOUR M. 


Data for the Blade 


BOGDONOFF* 


Langley Memorial Aeronautical Laboratory, N.A.C.A. 


SUMMARY 


Low-speed and high-speed two-dimensional cascade tunnels 
and a low-speed test blower were used to determine an experi- 
mental method of designing high-performance axial-flow com- 
pressor blades. In these apparatus, the important parameters 
involved in compressor blade design, turning angles or loading, 
design angle of attack, and qualitative critical speed have been 
evaluated. 


SYMBOLS 

Cl = blade lift coefficient referred to mean air 

M = Mach Number, ratio of velocity entering the cascade 
to the speed of sound 

M.,, = entering Mach Number at which sonic velocity is 
reached on the blade 

Mp = entering Mach Number at which losses through the 
cascade suddenly increase 

a = local pressure coefficient over the blade 

P., = pressure coefficient corresponding to sonic velocity 


for a particular entering Mach Number 
P;/P, = static pressure ratio across the cascade 


Q = quantity flow, cu.ft. per sec. 

U = rotational velocity of rotor blade element, ft. per sec. 
(Fig. 1) 

Va = axial velocity, ft. per sec. (Fig. 1) 

W = velocity of air relative to rotor element, ft. per sec. 
(Fig. 1) 

a = angle of attack, angle between effective entering air 
and blade chord line, deg. 

B = stagger angle, angie of effective entering air measured 
from axial direction, deg. (Fig. 1) 

n = adiabatic rotor efficiency 

6 = effective angle through which air is turned relative 
to the rotor, deg. (Fig. 1) 

o = solidity, blade chord divided by the gap between 


blades (Fig. 1) 
Subscripts 
0 = mean conditions, one-half vector sum of entering and 
leaving conditions 
1 = entering rotor 
2 = leaving rotor 


INTRODUCTION 


Tp RECENTLY, most aircraft gas turbines have 
been designed using centrifugal compressors, 
which are light, are easy to build, and give high pressure 
ratios per stage. In addition, design information was 
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already available from numerous supercharger tests. 
As the development of the gas turbine progressed, how- 
ever, the axial-flow compressor came into prominence 
because of the advantages of high efficiency and small 
outside diameter for large mass flows. The main dis- 
advantage, the low pressure rise per stage that has been 
obtained in commercial machines of this type (on the 
order of 1.18), necessitates the use of numerous stages 
to obtain the pressures that may be used efficiently in 
gas turbines. Many stages mean large numbers of 
accurately constructed blades and high weight. 

For the last three years, the N.A.C.A. has been work- 
ing on a program to develop axial-flow compressor blad- 
ing to give considerably higher pressure rises per stage 
than those now obtained in commercial machines and, 
at the same time, to give very high efficiencies. The 
theoretical approach had not yielded any practical re- 
sults that could have been applied directly to compres- 
sor blade design, so an experimental investigation was 
carried out by the Gas Dynamics Section of the Langley 
Memorial Aeronautical Laboratory of the N.A.C.A. 
This has been accomplished primarily by obtaining 
fundamental data on the maximum efficiency, loading, 
and operating Mach Number of a series of compressor 
blade sections. 

The program was divided into three parts: (1) A 
series of low-speed two-dimensional cascade tests were 
made to provide basic information about the turning 
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Simple vector diagram for a compressor rotor. 
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Fic. 2. Vertical cross section of two-dimensional low-speed 
cascade tunnel. 


angle and shape of the pressure distributions—the first 
to provide accurate data for the determination of power 
input and the second to show separation over the blade 
and provide a rational means for choosing an efficient 
operating condition; (2) low-speed rotating tests were 
made to obtain efficiency, to check the cascade tunnel 
turning angle data, and to find the maximum loading 
that could be used successfully in an actual rotating 
machine; and (3) these basic data were extended to 
high Mach Numbers to find effects of compressibility 
on turning angles, design angles of attack, shape of the 
pressure distributions, as well as to find the Mach Num- 
ber for critical speed and forcebreak. 

From the results of these three series of tests, a few 
typical sets of curves are presented which give the two- 
dimensional design data, the efficiency variation with 
loading, and the critical speed of the blade sections. 
As an example of the application of these results, an 
analysis was made to show the geometry—.e., the stag- 
ger angle and blade hub-to-tip ratio—for an optimum 
machine based on maximum pressure rise and on maxi- 
mum power absorption. 


Low-SPEED Two-DIMENSIONAL CASCADE TESTS 


The first of the investigations'~* was carried out in 
the low-speed two-dimensional cascade tunnel shown 
schematically in Fig. 2. The flexible walls at the top 
and bottom of the tunnel were so adjusted as to simu- 
late the flow through an infinite grid, although only five 
airfoils were used. Static pressures were measured 
ahead of the blades, and pressure distributions were 
‘taken over the center blade of the cascade. In ad- 
dition, yaw surveys were made behind the cascade to 
determine the turning angle. The tunnel was driven 
by a 25-hp. blower, which kept the region ahead of the 
settling chamber slightly above atmospheric pressure. 
The velocity through the test section was approximately 
100 ft. per sec. The blades tested were N.A.C.A. 65- 
series airfoils with the trailing edge slightly thickened 
to make a practical blower-blade section. The cambers 
varied from zero to that for a lift coefficient of 1.8 in 
free air. 
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A typical set of pressure distributions for a range of 
angles of attack are shown for one of the high-cambered 
blower blades in Fig. 3. The pressures over the airfoil 
are plotted as pressure coefficients, the local pressure 
minus the entrance pressure divided by the entrance 
dynamic pressure. Eastman N. Jacobs, who started 
this program, proposed, on the basis of isolated airfoil 
tests, that the minimum losses occur at the angle of 
attack at which the pressure distribution shows no 
peaks or has a uniform loading. This point is desig- 
nated the design point or design angle of attack—in the 
case shown in Fig. 3, it is approximately 17°. 

Tests of this type were conducted for a series of air- 
foils with entrance conditions corresponding to staggers 
of 45°, 52.5°, and 60° and solidities of 1.0and 1.5. The 
final results have been correlated in a set of compressor 
blade design charts, one of which is shown in Fig. 4. 

From the charts, for any turning,angle, stagger, and 
solidity, a blade section and corresponding angle of 
attack may be chosen which will give the desired per- 
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Fic. 3. Low-speed pressure distributions over N.A.C.A. 65- 
(12)10 blower blades. 8 =45°,o = 1.5. 
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Fic. 4. Axial-flow compressor blade design charts. 
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formance with a peak-free pressure distribution over the 
blade. For example, it is desired to turn air 20° at a 
stagger of 45° and a solidity of 1.0. The design camber 
is then about 1.1 (design lift coefficient of the isolated 
airfoil) and the angle of attack, 14°. 


Low-SPEED Rotor TESTS 


The second part of the investigation, the low-speed 
rotor tests,‘ was carried out in the low-speed single- 
stage test blower shown schematically in Fig. 5. The 
rotor was driven directly through a long extension shaft 
by a 75-hp. variable-speed motor. The air entered 
eight entrance ports, whose area could be varied by 
sliding ports, and passed through three screens, the 
converging entrance section, the test section, and then 
through an annular diffuser to the atmosphere. The 
test machine had a hub-to-tip ratio of approximately 
0.78. 

Since the main purpose of these tests was to show the 
applicability of the cascade results, only a single rotor 
was investigated. Complete surveys of total pressure, 
static pressure, and yaw angles were made at positions 
one-half chord ahead of, and behind, the rotor. A 
photograph of the actual test machine with the top half 
of the test section removed is shown in Fig. 6. The 
survey instruments and their positions are shown in 
Fig. 7. Each instrument consisted of a total-head tube, 
static-pressure tube, and yaw head on the same stem 
and provided for simultaneous readings of all three. 
The ‘‘null’’ method of taking data was used—that is, 
the tube was rotated until the yaw head showed the 
instrument was facing directly into the air stream. 
The total pressure, static pressure, and yaw angle were 
then noted. The tests were made with five sets of 
blades that were designed to set up free-vortex flow 
with lift coefficients at the pitch section varying from 
0.3 for blade 1 to 1.2 for blade 5. The solidity at the 
pitch section was 1.0, and the blade chord was held 
constant. 

Fig. 8 shows the variation of the turning angle along 
one of the blades. The points indicate the measured 
values; the dotted lines, the predictions from cascade 
tests. As can be seen, the agreement is good over the 
portion of the blade not affected by the hub and tip, the 
maximum deviation being less than 1°. 
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Fic. 5. Schematic diagram of single-stage test blower. 
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Single-stage test blower with top of test section re- 
moved. 


Fic. 7. Installation and close-up of survey instruments and 
measuring head. 
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Fic. 8. Comparison of measured and predicted turning angles. 
Blade No. 3. 


The efficiencies of the blades tested are shown in Fig. 
9, where the blade loading is represented by the lift 
coefficient at the pitch section. It is at once obvious 
that there is no decrease in maximum efficiency with 
increased loading. Instead, the peak efficiency rises 
sharply as the loading is increased and then remains 
constant with loadings up to the highest tested. The 
curve showi@g the variation of efficiency with quantity 
coefficient (the volume flow divided by the product of 


of 5 

The 
ssor 

A > " " 
) 4 
65- q 
— | 

29 


LIFT COEFFICIENT, 94 !9 
Clo 


1.0- 


0.6- 


= LIFT CURVES 


CONSTANT EFFICIENCY CONTOURS 


ANGLE .OF ATTACK 
Fic. 9. Variation of efficiency with loading. 


100- 
95- 
DESIGN POINT 
90- 
85- 
80 T T 
.50 .60 .70 .80 


QUANTITY COEFFICIENT 


Fic. 10. Variation of efficiency with quantity coefficient. 
Blade No. 3 


Va, 


Va,+ Vao 
2 


Vao 


Vector diagram for compressor rotor with varying 


Fic. 11. : 1 
axial velocity. 


rotor speed and tip diameter cubed) is presented in Fig. 
10. The condition for which the blade was designed is 
marked by a vertical line that is close to the peak 
efficiency. 

These tests show conclusively that the data obtained 
from the low-speed cascade tunnel can be used in the 
accurate design of low-speed compressor@Mlades loaded 
much higher than those now in use. The turning angle 
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checked the design turning angle within 1°, and the 
occurrence of peak efficiency at the design point verifies 
the assumptions that were made in the original cascade 
study. 

A further investigation showed that the two-dimen- 
sional cascade data could also be applied to blades 
incorporating radial flows if the blades were designed for 
equilibrium conditions from a modified vector diagram 
based on the average axial velocity. This modified 
diagram, on which the effective stagger and turning 
angles are defined, is shown in Fig. 11. Use of these 
effective angles to design blades for changes in axial 
velocity up to 40 per cent of the initial axial velocity 
gave agreement within 1°—2° between experimental and 
predicted turning angles. The efficiencies were as high 
as those for the two-dimensional or free-vortex blading. 


HIGH-SPEED Two-DIMENSIONAL CASCADE TESTS 


The third part of the Gas Dynamics program, that of 
extending the basic data to high Mach Numbers, was 
carried out in a high Mach Number two-dimensional 
cascade tunnel similar to that used in the original low- 
speed investigation. The tunnel is shown schemati- 
cally in Fig. 12. Because of the excessive power re- 
quired, no boundary-layer control was incorporated in 
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Fic. 12. Diagram of high Mach Number two-dimensional cas- 
cade tunnel. 
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this tunnel. High-pressure air from a large tank was 
exhausted through the cascade tunnel to the atmos- 
phere. Pressure distributions were obtained over the 
center blade, static pressures were measured ahead of 
the cascade, and the turning angles were measured by a 
yaw survey and photographs of tufts in the air stream. 
The tests were repeated with the center three blades 


mounted between glass plates so that schlieren photo- 


graphs could be obtained. 


A typical set of pressure distributions over a range of 
Mach Numbers is shown in Fig. 13, arranged in order 
of increasing Mach Number. As the Mach Number 
increases, the pressure distribution over the blade 
changes in a manner similar to that experienced with 
isolated airfoils. Peaks occurred on the nose of the 
airfoil, and these peaks increased until the critical speed 
was reached. If the Mach Number is increased still 
further, large shocks are noticed. 


In Fig. 14, schlieren photographs of the same blade 
section are shown, again arranged in order of increasing 
Mach Number. In the first pictures, the critical speed 
has not yet been reached, and no shock waves are notice- 
able. It should be noted here that what appears to be 
a thick blade boundary layer in these photographs is 
actually the flow at the corner of the blade wall inter- 
section. As the Mach Number is increased, small 
shocks occur and grow stronger and, finally, separation 
occurs (noticeable by the mottled area over the rear of 
the blade). This separation or change in the flow is not 
noticeable until some time after the critical speed has 
been exceeded. 


Fig. 15 shows the variation of turning angle with 
Mach Number for two blower blade sections. The 
upper curve shows a typical case; the lower curve shows 
the case for maximum change in turning angle before 
forcebreak, approximately 1°. At forcebreak a sharp 
drop occurs as the lift suddenly decreases because of the 
separation. 
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Fic. 15. Variation of turning angle with Mach Number. 
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Fic. 16. Variation of pressure ratio with Mach Number. N.A.- 
C.A. 65-(12)10 blower blades. a = 14.99; B =45°; o =1.5. 
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The curve of pressure ratio across the cascade versus 
entrance Mach Number is shown in Fig. 16, with the 
critical speed, M,,, and the inflection point in the curve, 
or forcebreak, My, designated. Since the large losses 
due to the boundary layer on the walls decrease the 
pressure ratio, this curve does not truly represent the 
pressure ratio that would be expected across a row of 
blower blades. However, any sudden increase in the 
losses above the boundary-layer losses on the walls 

. Should show up as a break in the curve. It is at once 
obvious that the critical speed is exceeded before force- 
break occurs. From these tests, it appears as though 
the critical speed is not indicative of either increased 
losses or changes in performance. Variation of the 
Mach Numbers for critical speed and forcebreak are 
shown versus the angle of attack for a typical blade 
section in Fig. 17. The design angle of attack of the 
section is marked by the vertical line. It can be seen 
that the peak Mach Number for critical speed and force- 
break is close to this design condition. Therefore, the 
low-speed design point is also the high-speed design 
point, and the low-speed turning angle may be used up 
to forcebreak. 


APPLICATION OF RESULTS 


The actual results of this series of tests apply only to 
the particular set of blades tested, but the apparatus 
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will, in addition, show the most efficient blade section if 
a series of tests of other blade sections is made. The 
data obtained are of primary importance in the estab- 
lishment of a set of limiting conditions for the maximum 
loading, efficiency, and Mach Number that can be ob- 
tained with any series of compressor blades. These 
limits establish the maximum pressure ratio or the maxi- 
mum power that can be obtained in a single stage of a 
compressor. A theoretical analysis of the variables 
involved in compressor design then becomes powerful 
and practical. For example, a two-dimensional analy- 
sis using the maximum turning angles and correspond- 
ing critical speeds from the low- and high-speed cascade 
tunnels shows that the maximum pressure rise will be 
obtained at high stagger angles, while the power input 
reaches a maximum at about 55°-60° stagger angle. 
If this analysis is extended to three-dimensional vortex 
flow compressors, the results show maximum pressure 
rise as the blade hub-to-tip ratio approaches 1.0—that 
is, short blades. The maximum power input, however, 
occurs when the blade hub-to-tip ratio is approximately 
0.7. 

For any specific application, the fundamental infor- 
mation on which to base the geometry of the compressor 
is now available. The limiting conditions will no doubt 
change as additional research is completed, but, with 
the loadings and Mach Numbers shown possible by 
these tests, performance comparable to the best existing 
axial-flow compressor could be obtained with half the 
number of stages. 


CONCLUDING REMARKS 


The results obtained in this investigation showed 
that: €1) the low-speed rotor and high-speed cascade 
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tests verify the assumptions that formed the basis for 
the design method based on low-speed cascade tests; 
(2) the turning angles, found from low-speed cascade 
tests, vary less than 1° for all Mach Numbers below the 
forcebreak; (3) maximum efficiency was obtained at 
the design conditions, and these conditions do not 
appear to change with Mach Number; (4) blades de- 
signed from these data have high efficiency and wide 
operating range, even for blades loaded much higher 
than those now in use; (5) the critical Mach Number 
must be definitely exceeded before the performance of 
a compressor blade changes radically, in a manner 
similar to the isolated airfoil at high Mach Num- 
bers. 

Tests of this type give the much-needed limiting con- 
ditions that establish the maximum pressure ratio or 
power input that can be obtained in a single stage of a 
compressor. With the loadings and Mach Numbers 
shown possible by these tests, performance comparable 
to the best existing axial-flow compressor could be . 
tained with half the number of stages. 
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Letter to the Editor 


Dear Sir: 


In the paper ‘‘Computation of Influence Coefficients for Air- 
craft Structures with Discontinuities and Sweepback” by Samuel 
Levy (JOURNAL OF THE AERONAUTICAL ScIENCES, Vol. 14, No. 
10, p. 547, October, 1947), a method of representing the distrib- 
uted mass of a wing by a system of concentrated masses is pre- 
sented. A six-mass approximation is given. The location of 
these masses can be arbitrarily chosen. Their magnitudes are 
determined from a system of six simultaneous equations. It has 
been noted that, when these six masses are located on two 
straight lines, the equations become incompatible, and a formal 
solution is impossible. In order to eliminate this difficulty, the 
following procedures are suggested. 


A nine-mass approximation will be treated. If we have an 
airplane wing in mind, it would be convenient to locate these 
masses on the ribs and the beams of the primary structure, since 


the influence coefficients for these points can be readily evaluated. , 


The system of points is shown in Fig. 1. The X-axis is taken to 
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be parallel to the fuselage axis, and the Y-axis is directed toward 
the wing tip. The ribs are assumed to be parallel to the fuselage 
axis, which is true in most cases. The wing will be divided into a 
finite number. of strips. One of these strips is shown in Fig. 1. 
In fact, the origin of the coordinate axes shown coincides with 
the centroid of the strip. Let the maximum displacement of 
J this strip be approximated by the second-degree surfaces. By 
‘means of similar reasoning as in Mr. Levy’s paper, it is possible 
to specify the following conditions to be satisfied by the nine 
concentrated masses: 
= S dm 
= 
= JS dm x? 


= 
2 Mijxizyiz? = S dm xy? (i,j 1, 2, 3) (1) 


2M = S dm x* 
= S dm 
= 
Since the y coordinates are the same for the points P,;, P2;, 
and P3;, it is possible to simplify the solution of Eq. (1) by the 
following scheme. Introduce the expressions 


M; 1 1 1 Mi; 
M;%; = My (2) 


Eq. (1) can be reduced to the following: 


Mz Mopy? = M (3) 
ys? Ms Moapys? 


where 

fdmx Sfdmx* 

Sdmy Sfdmxy 

Sdmy? Sdmxy? fdmx*y 

Denote the inverse of the first matrix in Eq. (3) by A, the solu- 
tion can be written as 


M = 


1 M, M, Pyi? 
Moy? | = (4) 
M. Mk; Mspys? 
where 
An 
A = Ay n Ya 
Ay 
From equation (2) it follows that 
Mi; M; 
M3; = Asi (Gj = 1, 2, 3) (5) 
M3; Mipyi* 
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where 
1 1 1 ithe 
4:3 = X2j 
az? 


The magnitude of the concentrated masses is therefore given by 


Mi; 
M2; = 
M3; 


where the prime denotes the transposed matrix. 

If matrix iteration method is used to evaluate the vibration 
modes, the mass matrix will have only diagonal elements which, 
in turn, contain contributions from all strips on the wing. 

A better approximation can be obtained by the consideration 
of the energy involved in free vibration. The kinetic energy of a 
continuous system can be approximated by a system of general- 
ized masses associated with an arbitrary set of points. In the 
case of nine points, the surface through these points can be repre- 
sented approximately by 

z= NA= 

a+ by + cy? + dx + exy + fxy? + gx? + hx*®y +hx*y? 


(j = 1, 2, 3) (6) 


where 
A = {a,b,c,d,e,f. g, h, k} 
N = y, y?, x, xy, xy, 2%, x%y, x%y?} 
are two column matrices. The kinetic energy involved will be 
T = fdm(2/2) = A'MA/2 (7) 
where M = {dm NN’. 


The coefficients a, b, etc., can be expressed in terms of the maxi- 
mum displacements at the nine arbitrary points. Thus, 


A= AZ (8) 
where Z = {21, 22, 23, ... 2} is a column matrix of displacements 


at the nine arbitrary points. Substituting Eq. (8) into Eq. (7), 
the kinetic energy is reduced to 


T = 
and the generalized forces acting at the arbitrary points are 
Q = (d/dt)(dT/dZ) = [A’MalZ (9) 


The matrix product [A’M/A] is the generalized mass matrix as- 
sociated with the arbitrary set of points. It is different from the 
diagonal matrix discussed previously in that coupling terms do 
appear. Therefore a better approximation to the actual dy- 
namic system can be expected. 


Conrap C. WAN 
Graduate School 
Illinois Institute of Technology 
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Range Performance of Turbojet Airplanes 


IRVING L. ASHKENAS* 
Northrop Aircraft, Inc. 


SUMMARY 


Expressions are developed for the maximum range of a turbo- 
jet-propelled airplane for three cruise flight conditions: constant 
altitude, varying altitude below 35,000 ft., and varying altitude 
above 35,000 ft. It is shown that maximum range is obtained by 
increasing the cruising altitude as the fuel is consumed and that 
this range is approximately given by the derived expressions even 
under conditions of nonoptimum operation. It is further shown 
that the thrust loading, 7/W, is of prime importance, the range 
for a given airplane increasing with this parameter at a rate some- 
where between the square root and the first power of this ratio. 

Airplane configurations are derived which correspond to the 
aerodynamic optima for the cruise conditions investigated. It is 
indicated that, except for extreme cases, these configurations may 
be unreliable guides to optimum airplane designs. 


INTRODUCTION 


§ be RANGE EQUATIONS for turbojet-propelled air- 
planes differ basically from those of conventional- 
engine aircraft in that they are based upon a fuel con- 
sumption that is approximately proportional to the 
engine thrust rather than to the power. The basic 
equation is thus written: 


-—dR = VdW/cT 


where R is the range in miles, W is the airplane weight in 
pounds, T is the thrust in pounds, c is the specific fuel 
consumption in pounds per hour per pound of thrust, 
and V is the speed in miles per hour. Assuming the 
normal analytical expression for the drag polar and 
equating thrust t8 drag (level flight), 


—(dR/dW) = (2/cpSV)[Cop + 


S being the airplane wing area, p the atmospheric den- 
sity, Cp, the airplane parasite drag coefficient, C; the 
lift coefficient, and A, the “effective’’ aspect ratio.! 
Substituting for the velocity in terms of C;, 


Re (oy 
1 


RANGE AT A CONSTANT ALTITUDE 


If now we assume a constant altitude, then p becomes 
constant and c may be assumed essentially constant so 
that the maximum value of dR/dW (miles per pound) is 
obtained at the minimum value of the expression in 
parentheses. It is easily verified that this occurs when 


C,?/ TA, f (2) 


and the maximum value of dR/dW is 
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dW 4c YpSW\ 3 
The formula for maximum range is obtained by inte- 


grating from the initial weight, Wi, to the final weight, 
W2, and, with the proper constants, can be written as 


_ 30 ( 
= 1 W, (miles) (3) 


where o is the atmospheric density ratio. 

From the above expression it is clear that the range 
will steadily increase with increasing altitude, so that 
the maximum economy for any given weight will occur 
at an altitude corresponding roughly to the “cruise 
ceiling’”’ for that weight. The maximum overall range 
will thus be obtained by climbing steadily as the fuel 
weight is consumed. The rate of climb required to 
maintain this optimum altitude-weight relationship 
will generally be so small as to be negligible; it is there- 
fore possible, for this cruise climb, to equate thrust to 
drag with little error. é 


RANGE AT VARYING ALTITUDE 


If we do this and assume that the thrust variation 
with altitude can be represented as some power, m > 0, 
of o, 

Cp ) 
Too” = — W WwW 
us CL + TA, 
where 7) represents the sea-level value of the cruise 
thrust being used. This equation defines the altitude, 
attainable at a given rating, as a function of weight. 
By substituting it in the basic Eq. (1), we can write, 


dw ic poSW W CL 4A, 


This expression is intended to show that, for a given air- 
plane and an approximately constant specific fuel con- 
sumption, the value of dR/dW obtainable in a climbing 
cruise increases as the thrust setting is increased. This 
result is self-evident from the basic consideration that 
the “‘cruise ceiling’’ will increase with thrust rating. 


Engine Performance Versus Altitude 


Now, the rated (cruise or normal) power-plant char- 
acteristics may be approximated by: for altitudes be- 
low 35,000 ft., 
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for altitudes above 35,000 ft., 
T = T, 35(0/ 035) 


C = 


where the subscripts refer to sea-level and 35,000-ft. 
altitude conditions. The above approximations have 
been checked against a number of American and British 
turbojets with reasonably good results, provided the 
equations are applied to velocities in the region from 300 
to 600 m.p.h. 

The maximum range, which corresponds to a “‘cruise- 


ceiling” operation at,the maximum allowable continuous 


power as discussed above, may readily be calculated 
using the above expressions. 
Maximum Range Below 35,000 Ft. 


We may write then, for the portion of the range which 
is flown below 35,000 ft., 


dR 1 


Co, 
155 + 
dW co? SW Cr TA 


This equation is clearly obtained by substituting for c 
in Eq. (1). Now as before, we determine the cruising 
altitude as a function of the weight by equating thrust 
to drag. Thus, 


= [(Cp,/Cr) + |W (4) 


and, substituting for 0°’, the range equation becomes 


For maximum range, minimizing the expression in 
parentheses yields 


= */sCoy (6) 
and the maximum value of dR/dW becomes 


aw 64\5 poSW? 


1.00 
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Fic. 1. Range reduction for nonoptimum operation. 


Integration and conversion to the proper constants 
gives 


R= 5/4 val - 


Maximum Range Above 35,000 Ft. 


In an exactly analogous manner we can calculate the 
maximum range attainable above 35,000 ft. Thus, 


= + (8) 
and, by substituting this stm for o in Eq. (1), 


aR 1 27% Cp, 
dw C35 o35P09 CL TA, 
Minimizing the expression in parentheses gives 


= (10) 


Substituting this result in Eq. (9) gives, for maximum 


dR/dW, 
dw \3 pS Cp, 


and, finally, 


24.2 (*) 
= S Cy TC, W, (miles) (11) 


Maximum Total Range 


1) (miles) (7) 


To calculate the maximum range exclusive of the 
initial climb and final descent, it is necessary to deter- 
mine the weight corresponding to the optimum cruise at 
35,000 ft. This figure becomes W2 in Eq. (7) and W; in 
Eq. (11), the initial and final airplane weights being W; 
and W2, respectively, in these equations. There is some 
difficulty in determining the 35,000-ft. weight, since 
Eqs. (4) and (8), when the corresponding optimum value 
of C; is substituted in each, yield values for the weight 
which are different by about 3 per cent. This effect is 
negligible on the total calculation. 


EFFECT OF DEVIATIONS FROM THE OPTIMUM CRUISE 
SPEED 


It appears, in general, that the maximum range is 
obtained for the three cases noted above by flying at a 
constant C;, corresponding to C,?/7A, = kCp,, where k 
varies with the flight conditions assumed. The depend- 
ence of the maximum range on the choice of k can be 
determined by examining the effect of variations in k 
on the constants appearing in the final range equations. 
It is easily verified that these are proportional to: 


+ constant altitude 
(k~*/* + k*)-? varying altitude below 35,000 ft. 
+ varying altitude above 35,000 ft. 


The effect of departing from the optimum, fo, on the 
range constants is shown in Fig. 1. It is clear that for 
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RANGE PERFQRMANCE OF TURBOJET AIRPLANES 


large variations in k the effect on range is reasonably 
small, being of the same order, for example, as in the 
case of the Breguet? equation for reciprocating-engine 
airplanes, for which the range may be shown to be pro- 
portional to (k~ 4 This means that the fore- 


going range equations may be used with good accuracy — 


even though the airplane cruising speed does not corre- 
spond to the optimum. Thus, for example, should the 
optimum cruising speed, calculated with no regard for 
compressibility effects, exceed the airplane critical 
Mach Number, the cruising speed might be reduced 
and the effect on the calculated maximum range would 
nevertheless be small. This last conclusion contradicts 
the result of reference 3. 


IMPLICATIONS OF THE RANGE EQUATIONS 


The three expressions for range can best be compared 
by noting that (Vx/ 2)V A./Cp, is equal to the 
maximum lift-drag ratio (L/D)maz. and by defining 


= (To/W,i)(L/D) maz. 
Then Eqs. (3), (7), and (11) may be written as: 


_ 31.9 ( 


‘ea’, ( (7a) 


Co SCp, \D W: 
_ 23.0 Wi (Z)" 
R= n SCop \D) mar log, (11a) 


The airplane parameters, other than weight ratios, in- 
volved in these equations may be expressed as: 


V Ws/SCo,(L/D) maz." 


where x varies from zero to one to one-half as the cruise 
conditions are changed. It is now again apparent that 
the thrust loading, 7)/W1, is a parameter of prime im- 
portance for long-range aircraft and that, in fact, the 
range belew 35,000 ft. for a given initial weight varies 
directly with the first power of this parameter and above 
35,000 ft. with the square root of this parameter. This 
analysis emphasizes the fact that in expressions (3) the 
value of o is determined, essentially, by the value of 
T/W. 

For airplanes with normally high thrust loadings, the 
entire flight can be made above 35,000 ft., and the opti- 
mum true speed will then be constant. Expressing the 
range in terms of this optimum cruising velocity, 
Eq. (11) becomes 


R= (0.94/c3s5) Vo(L, log, We) (miles) (11b) 


where Vo is the cruising speed in miles per hour. This 
expression demonstrates that, for a high-performance 
turbojet airplane, increasing the design cruising speed 
will reflect a similar increase in maximum range. High 
speed and range are thus compatible and, indeed, com- 


IAN 


Fic. 2. Effect of me area on optimum range. 
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plementary requirements, and wing loading is not, as 
indicated in Eq. (3), a prime parameter. 


OPTIMUM AERODYNAMIC CONFIGURATIONS 


While wing loading, as such, has a secondary effect 
on range, wing size (and therefore wing loading) will be 
reflected in the magnitude of Cp,, which influences the 
range as shown in the various formulas. It is thus 
possible to calculate an optimum wing loading or, more 
correctly, an optimum relationship between wing and 
fuselage drag. 


Given Fuselage-Nacelle Parasite Area 


Suppose that the parasite area is given by a constant 
part fi, necessary for a given power plant, and the wing 
parasite area Cp,S—.e., 


Cop = Cor + (fi/S) 


From Eq. (3), best range at constant altitude for a given 
airplane (i.e., const. W, T, A,, etc.) occurs at the mini- 
mum value of 


[Coo + 


Differentiating with respect to S results in an optimum 
when 


S = 1/2( fi/Coo) (12) 


For range at varying (optimum) altitude below 35,000 
ft., the best range according to Eq. (7) will occur at the 
minimum value of 


[Coo + 
which corresponds to: 


S = */2(fi/Co») (13) 
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Fic. 3. The variation of maximum range with volume disposi- 
tion. 


and, finally, for range at varying altitude above 35,000 
ft., minimizing 


[Cos + (fi/S)]S'* 
[see Eq. (11)] gives: 
S = fi/Coo (14) 
In general, the wing area may be expressed as 
S = n(fi/ Coy) 


The pertinent functions of m defining its effect on the 
range are: 


+ constant altitude 
+ varying altitude below 35,000 ft. 
(n'/* + n~)-1 varying altitude above 35,000 ft. 


These functions are shown in Fig. 2, which indicates 
that fairly large variations from the optimum wing area 
will affect the range but slightly. Under these circum- 
stances and considering the fact, as previously men- 
tioned, that, for high-performance jet airplanes, the 
best range will occur above 35,000 ft., Eq. (14), S = 
fi/Co, Tepresents a good compromise for an aero- 
dynamic optimum, under the conditions outlined above. 


Given Volume of Airplane 


Suppose now that the airplane is required to enclose a 
certain volume, v, which may be distributed betweef 
the wing, v,., and fuselage, v,. In that case the mini- 
mum drag coefficient may be written as: 


Cop = C2 + 


where C; is the effective skin-friction coefficient and 
where S,, the fuselage wetted area, is equal to k,v,”*. 
Expressing 


S = 


‘and substituting for the functions of Cp, and S, which 


correspond to the best range (see above), the range may 
be shown to be proportional to: 


{ + (Rn/ Rw) 


varying altitude above 35,000 ft. 

The variation of this function with v/v, and k,/k, 
from the value at v/v, = 1 is shown in Fig. 3. It is 
clear that the optimum v/v,, giving a minimum for the 
ordinate, depends on the value of k,/k, being consid- 
ered. Thus, for values of k,,/k, less than approximately 
0.275, a wing-fuselage configuration, with the relative 
fuselage size increasing as k,/k,. decreases, is indicated. 
On the other hand, for values of k,/k,. greater than 0.275, 
the all-wing configuration gives best results. 

For bodies of revolution, k, may be approximated by 
considering two extreme shapes: a cylinder with hemi- 
spherical ends and a double cone. For the former, 
ky, = 3.68 (i/d)", approximately, and for the latter, 
3.84 (I/d), I/d being the fuselage length-diameter 
ratio. For wings of elliptic section : 

Ry = 2.25[A'*/(t/c)”"] 


approximately, and for wings of double-wedge section 
304 [A'/*/(t/c)"], A being the geometric aspect ratio 
and (t/c) the airfoil thickness ratio. A reasonably 
accurate value for k,,/k,. may be obtained by averaging 
the above figures, whereby 


kn/Rw = 1.44 [(t/d)/A 


= 
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Fic. 4. Effect of airplane configuration on value of kn/kw. 


I 
Cc 
- p 
e 
ir 
pat 
lo 
1: 
0] 
be 
fo 
Tl 
va 
tre 
fig 
OV 
- ‘a 
ver 
dy 
tail 
I 
TICs 
erro 
4 -50 plat 
T 
wert 
text 
— 
orig’ 
lette 
5 
Dr. 
= jf 
Inst 
Fe 
| siste 
ther 
—— 
bein 


and 


lich 
nay 


/kew. 


RANGE PERFORMANCE OF TURBOJET AIRPLANES 101 


This function is plotted in Fig. 4 for two values of the 
ratio, (1/d)/A. The configurations sketchedarebasedon 
current structural practice and correspond to the maxi- 
mum wing aspect ratio attainable for the thickness ratio 
at which each is plotted. Inspection of this figure 
shows that for any wing greater than 12 per cent in 
thickness, the value of k,,/k, for suitable aerodynamic 
configurations will be greater than 0.275. This value 
makes no allowance for structure, control surfaces, 
power plant, etc., but assumes that the entire volume of 
both wing and fuselage are available for pay load and 
equipment. Under these assumptions it appears that, 
in general, for wings greater than 12 per cent in thick- 
ness, the all-wing version will be the ‘‘optimum enve- 
lope’ aerodynamically, whereas for wings thinner than 
12 per cent, the wing-body configuration will be an 
optimum. 


CONCLUDING REMARKS 


The foregoing optima have been calculated on the 
basis of aerodynamic requirements only and are there- 
fore to be considered with a great deal of suspicion. 
This is especially true since it has been shown that 
variations from the optimum produce, in general, ex- 
tremely small changes in the range; the effect of a con- 
figuration change on the weight ratio could thus easily 
overpower any beneficial aerodynamic effect.* The 
“optimum envelope’ calculation comes closest to elimi- 

* Thus, for a propeller-driven airplane, where range is in- 


versely proportional to Co,'/ *, the all-wing is the obvious aero- 
dynamic optimum but only becomes a true optimum under cer- 


tain conditions of size and loading. 


nating this objection, since it does take account, at least, 
of the volume available for pay load. For the extreme 
values of k,,/k,. shown on Fig. 3, for which the variation 
of range with volume disposition is large, it appears 
likely then that the aerodynamic optimum may come 
close to defining the best airplane for range. This 
means that thin-winged airplanes will be of the wing- 
body type and thick-winged airplanes will be of the all- 
wing type. 

The above cautionary note does not, of course, extend 
to the basic range equations herein developed. They 
are accurate expressions of the best range available for 
the power-plant and drag characteristics assumed. 
The fact that the incompressible drag equation has 
been used does not critically limit the analysis because, 
as has been shown, the equations will be accurate within 
5 per cent provided the optimum cruising speed does 
not exceed the critical Mach Number for drag diver-, 
gence by more than 15 per cent. In the case of air- 
planes that are severely Mach-limited, many of the 
above considerations will, of course, be invalid. For 
example, the form of Eq. (11b) shows that, for such an 
airplane, increasing 7/W will not, as indicated in 
Eq. (lla), increase the range. 
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Letters to 


Dear Sir: 
In a recent Letter to the Editor (JoURNAL OF THE AERONAU- 
TICAL SCIENCES, Vol. 14, No. 11, pp. 625-626, November, 1947), 


errors were made in discussing the subject of length effect, flat . 


plates on elastic supports, compressed. 

The methods of A. J. Miles, cases 1 and 2 of his original paper, 
were reproduced in Theory of Elastic Stability, 1st Ed., by S. 
Timoshenko. In so doing, notational changes were made in the 
text (correctly, as I have found), but the diagrams of the be- 
havior of y versus ¢? as published apply to the notation of the 
original paper. For the text of the treatise, the value ¢* should 
be changed to @. This error was not noted in my original 
letter, so the comments therein are in error, even though it 
appeared at the time that some 69 column tests carried out by 
Dr. George Winter at Cornell for the American Iron and Steel 
Institute were satisfactorily consistent when so reduced. 

Fortunately, the correct procedure is simpler, and the con- 
sistency of the interpretation of the column tests is improved 
thereby. Unfortunately, the tests now show little effect of the 
sort expected by Miles because of the shortness of the columns 
that lay in the range (L/r) between 23 and 62, most of them 
being below about 35. 


the Editor 


The columns comprised two equal channels back to back, 
usually bolted, with no attempt at securing common action of 
the two channel web plates. Each channel was a symmetrical 
section with unstiffened lip flanges. 

Let the clear width of web plate be b, the thickness 4, and the 
length a. Find the area thereof, Ap = bh, and the moment of 
inertia thereof, Jp, taken about its center-of-gravity axis parallel 
to the long dimension 6. Find also the cross-sectional area of 
the channel, 2A 7, and the moment of inertia thereof, 2/7, about 
its own center-of-gravity axis parallel to the plate dimension } 
(normally the Y—Y axis). 

Similarly, let the clear width of the lips be w and their thick- 
ness h. 

According to Miles, so long as (a/b) does not exceed a value 
(a/b) 1, the effective width of the web plate for the compressive 
stress f, is that of a plate supported rigidly on all edges or (b/h), = 
10,327/Vf., while the effective width of a flange is (w/h). = 
3,672/ Vie. If the overall web width is D, the overall lip width 
is B, the effective area of one channel, A,, is 


A — ([(6/h) — (6/h).] + 21 ((B — 1)/h] — (w/h).})h? 
and the calculated ultimate load on a short column is Py. = 


/ Rw 
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fy2A, for the pair of channels where f, is the yield point of the 


material. 

In the region (a/b) is less than (a/b)1, Y = 2x. When (a/b) 
is greater than (a/b)z, m = 1,6, = w/(a/b)1, and, for values 
of @ which are stall, as 0.5 to 1.0 or less, y = +/26 approxi- 
mately. ” If E = 29,500,000 for steel, » = 0.30, 


0.91([r/Ip) = + (2Ar/Ap)] 


so that 


(a/b), = 0.6745V (Ir/Ip)/[1 + (2Ar/Ap)] 


When (a/b) is greater than (a/b)z, the effective width for the 
region (a/b) is lessthan (a/b), must be multiplied by (a/b), + 
(a/b) so far as the web plate is concerned. 

When (b/h) is greater than (b/h)., (b/h). is used, but, when 
less, (b/h) is used; similarly for (w/h). 

The above methods were applied to all 25 cross sections from 
which 69 test sections were derived. In only four cases did the 
length effect result in a reduction in effective web width, and in 
no case did the reduction in width exceed 24. In several cases 
there was a material width reduction due to (b/h), being less 


than (b/h). In many cases (at least half) the length effect 


reduced (b/h), but not to the extent of making it less than (b/h). 

The ratios (Py/Pyc) of the actual ultimate strength, P,, to 
Prue, the calculated strength, were calculated, and an average 
was found for each of the 25 sections. These ratios were then 
plotted against (L/r) for the section. 

In general, an average line through the points in (L/r) and 
(Py/Puc) (25, 1.00), (50, 0.700) represented the data fairly well. 
All values of (L/r) less than 25 gave ratios greater than 1.000. 
In the region with (L/r) greater than 50, the test results lay above 
the average, in spite of the fact that virtually all columns were 
reported to show low strengths because of lack of symmetry in 
the sections. The sections subject to the Parr and Beakley 
torsional failure criterion, (B/D) is greater than '/,B?/Lh, all 
showed values on the low side. The exceptions to the above 
rules were all of 12 U.S. Std. material having an average yield 
point of 34,500 Ibs. per sq.in. but a range in coupon test values 
of 7,500 Ibs. per sq.in., suggesting a greater than average strength 
of the material. 

It was believed that the plot of the allowable compressive 
stress to the allowable extreme fiber stress ratio would be of 
interest, the allowable extreme fiber stress f, being taken as 
f,/1.85, while the allowable stress 


. (P/A.) = 0.464f, — 


when (L/r) is less than 24,000/Vf,, 
or (P/A,) = 134,000,000/(L/r)? 


when (L/r) is greater than the limiting value. This curve gives 
0.850 when (L/r) = 0. In most cases the formula curve lay 
below the plotted points. Hence, it is felt that the methods of 
analysis are fully consistent with the tests of Dr. George Winter. 

It seemed desirable to calculate the value of (Z/r) for each 
of the 25 typical cross sections at which (a/b) would be equal to 
(a/b)x. The range was from 15.7 to 61.8. It happened that 
the ratios (b/h) were tabulated in an adjacent column of the 
table. It was noted that these columns, made of two channels 
with unstiffened flanges placed back to back, had the relation 
(L/r)t = (6/h), which was correct to within three units in the 
worst case and much closer in the average case. 

If the above were true, then a column with (L/r) = 120 
would have (b/h) = 120 if length effect is to be avoided. But, 
for a material with f, = 33,000 or f, = 18,000 Ibs. per sq.in., 
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the ratio (b/h). = 56, approximately. Hence, only 56/120 of 
the web are effective if the length effect is avoided, But suppose 
(b/h) = 56, approximately, to avoid reduction due to (b/h) 
values. Then such a column would have (a/b), = 56, approxi- 
mately; hence, the length effect will make only 56/120 of the 
web effective, a result the same as before. 

Obviously, when long columns are involved, stiffer flanges 
are demanded than can be had with unstiffened legs. 

In passing, it may be noted that Case 2 of Miles shows no length 
effect. Flanges of beams with a stiffener on one edge and deep 
web on the other have an effective width (w/h). = 11,990/Vf,, 
provided the stiffener is adequate. At (w/h) = (w/h)., the value 
of (d/h) for overall lip depth to thickness should be at least 
0.28 (w/h).. A greater depth is required when (w/h) is greater 
than (w/h),, but not so great as d = 0.28w. When (w/h) = 
3,672/ Vin, no lip is required. Hence, a straight-line relation 
of (d/h) to (w/h) exists between the two effective width values. 

For duralumin sheets, multiply the coefficients in the effective 
width formulas by 0.5822 = VV 10,000,000/29,500,000. 

It is possible that, in the region of reduced web width, the 
values of Jy and Jp, together with Ar and A p, should be reduced. 
This would lead to solutions by trial. Lacking adequate test 
data in these cases, no attempt has beer’ made to make such 
corrections. 

The column tests herein reported are now consistent with such 
beam tests as were made at Cornell, in view of the procedure 
laid down. No beams had long enough compression flat ele- 


‘ments to show a length effect. In other words, given some 


guidance, the theory of elastic stability has proved of great 
assistance in securing calculated values for ultimate strengths 
in the case of the use of thin metal sections. 

Since most columns would show length effects even in the 
range of (L/r) values appropriate to primary compression 
members, this length effect cannot be ignored in designing struc- 
tures involving sections of thin metal. , 

Epwarp ApAMs RICHARDSON 
Publications 
Bethlehem Steel Company 


Dear Sir: 

Doubtless you have had the three typographical errors in 
Chang’s Letter to the Editor (August, 1947, JouRNAL OF THE 
AERONAUTICAL SCIENCES, p. 456) pointed out previously. They 
concern Busemann’s series for supersonic flow. 

The first two, the use of PV? in Eqs. (4) and (5) instead of 
pV? are of small import, but the expression for C, could cause 
serious concern to anyone employing it. It should read 


A correction is also in order for the expression for Busemann’s 
series as presented in ‘“‘Aerodynamic Characteristics of Rec- 
tangular Wings at Supersonic Speeds,”’ by Bonney in the Febru- 
ary, 1947, JouRNAL, p. 110. The second term in Eq. (1) should 
read 20 csc 2yo, not 2y csc? 

ALAN PoPE 
Associate Professor 
Daniel Guggenheim 
School of Aeronautics 
Georgia School of Technology 
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The Disturbed Flapping Motion of 
Helicopter Rotor Blades 


HENRY PARKUS* 


Vienna Technical University 


SUMMARY 


The article deals with the flapping motion of a helicopter blade 
when disturbed by external forces (gust, variation in cyclic pitch). 
The investigation is based on a rectangular and untwisted blade 
with the flapping hinge on the rotor axis. The disturbance gives 
rise to nonperiodic oscillations, which are described by a linear 
differential equation with periodic coefficients. In solving this 
equation, a criterion is obtained which permits rapid determina- 
tion of the conditions of stability. Numerical evaluation shows 
that the disturbed motion is stable, though at a certain critical 
advance ratio partial destabilization sets in which increases with 
increasing flying speed. The influence of coupling between 
flapping angle and blade pitch is also considered. Practically, 
the disturbed motion is finished after one revolution of the rotor. 


SYMBOLS 
= flapping angle of the blade 
v = azimuth of the blade 
Mr = thrust moment of the blade about the flap- 
ping hinge 
Mw = moment of the blade weight about the flap- 
ping hinge 
J = moment of inertia of the blade with respect 


to the flapping hinge 
dy/di = angular velocity of the rotor 


@ = 

R = blade tip radius 

r = radial distance of blade element 

c = blade chord length 

B<1 = tip loss factor 

Ce = lift-curve slope of the airfoil 

p = air density 

7 = blade angle 

x = r/R 

m = V/wR = advance ratio (_V = flying speed) 
ny = inflow ratio 

k = factor of coupling between flapping angle 


and blade angle. k < 0 indicates a re- 
duction of the blade angle with increasing 
flapping angle 

Bp = forced flapping angle 

By = free flapping angle 

y = R*cCy’p/J= blade mass coefficient 


INTRODUCTION 


WwW" A HELICOPTER ROTOR BLADE is disturbed by 
a gust or by the altering of the cyclic pitch posi- 
tion, free oscillations of the blade are superimposed on 
the steady flapping motion about the hinge. These 
free oscillations are treated in this paper.’ 

The investigation is based on arectangular, untwisted, 


Received January 31, 1947. Revised and resubmitted July 23, 
1947. 


* Wissenschaftlicher Assistant an der Lehrkanzel f. Festig- . 


keitslehre (Prof. Dr. Girkmann). 


and rigid blade, with its flapping hinge on the rotor 
axis. A variation of the flapping angle with the blade 
pitch is also considered. 

It is assumed that angular velocity of the rotor and 
the flying speed do not vary in magnitude and direction 
during the oscillations. 


THE EQUATION OF MOTION 


The flapping motion of a rotor blade, with centrally 
located flapping hinges, is governed by the following 
differential equation” * 


(@°B/dy*) + 8 = (Mr — My)/Jw* (1) 
Using the notation 
Uy =x+yusiny 


up = \ — cosy — x(dB/dy) 
= + kB 


and making certain simplifications,” * one can write 
the thrust moment in the following form: 


Mr = R'wcC,'(p/2) (up + Bur)urxdx (2) 
The angle 8 may be divided into two parts, 
B = By + By (3) 


Here, 8, the periodic solution of Eq. (1), corresponds to 
the “‘forced”’ flapping motion and is represented by the 
particular integral of the inhomogeneous differential 
equation; while §,, the nonperiodic solution, corre- 
sponds to the ‘‘transient”’ flapping motion and is repre- 
sented by the general solution of the homogeneous 
differential equation. 
Substitution of Eqs. (3) and (2) in Eq. (1) with the 
abbreviation 
= ¥ (4) 


leads to the homogeneous differential equation 
for 8, (accents mean differentiation with respect 


toy): 
Bs 
+ sin + By + (uy cos - 


B? B? 
uky sin + wky cos 2y + sin 2y+ 
Bs B? 
yk yk 0 (5) 


The coefficients of this differential equation are periodic 
functions with the period 27. 
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The theory of these differential equations has been 
developed by Floquet.‘ The question of stability of 
the solutions, for certain restrictions on the coefficients, 
has been treated by Calamai.* 

Before carrying out the determination of the solu- 
tions of Eq. (5), it will be useful to give a short outline 
of Floquet’s theory.® It will be convenient to write 
Eq. (5) in the general form 


By" + By + = 0 (6) 


where p and g are periodic functions with the period 
2x. Placing p = O, one obtains the well-known differ- 
ential equation of Hill. 


FLOQUET’S THEORY 


Let the functions 6,(y) and f2(y) be a fundamental 
system of solutions of Eq. (6). Then, because of the 


periodicity of the coefficients, the functions + 


and $.(y + 27) are also solutions and can therefore be 
expressed linearly in terms of the fundamental system: 


+ = + (7) 


Floquet's theorem now states that there always exists a 
solution B(y) of Eq. (6) so that 


By + 2m) = sB (8) 


s being a real or complex constant. A(y) can also be 
expressed in the fundamental system: 


= + (9) 
From Eqs. (8) and (7) it now follows that the two con- 
stants Ci, C2 are related by the equations: ° 
Ci(au — s) + Cran = of (10) 
+ — s) = 0 


The equations have a nontrivial solution only when 


aun ay 


2 = 
5 (an + 22) + an dn 0 (11) 
Usually s is written in the form 
s= (12) 


v is called the characteristic exponent of Eq. (6). The 
two (in general distinct) roots s; and sz of Eq. (11) cor- 
respond to the two characteristic exponents »; and », 
respectively. 

When 5s; and sz are thus determined, the ratio of the 
two constants C; and C, follows from Eq. (10). The 
functions B; and £, which satisfy condition (8) are then 
calculated by means of Eq.°(9). 

Floquet has shown that these two functions, B, and 
B2, also form a fundamental system. Therefore, the 
general solution of Eq. (6) can be written in the form 


where ¢;(y) and ¢g2(W) are periodic functions with the 
period 

If function 8; is bounded for arbitrary initial condi- 
tions, then a stable solution is obtained, and the con- 
dition 

R(v) < Oor |s| < 1 (14) 


is fulfilled. Here, |s| denotes the modulus of the real 
or complex number s, and R(v) means the real part of v. 

Eq. (11) is invariant with respect to the particular 
fundamental system that is chosen. Therefore, a sys- 
tem defined by the initial conditions 


=1, = 0 
= 0, = 1 
may also be selected. Eqs. (7) now give, for y = 0, 
= au, = da 
= ay, = ax 


Hence, Eq. (11) becomes 
s? — + + W(2x) =0 (16) 
W is the Wrofskian of the fundamental system ;, fo. 
By a well-known formula of the theory of determinants 
W(2r) = W(0) exp. [— p(y) dy] 
Noting that W(0) = 1, Eq. (16) can be written in the 
form 


s? — + B2’(2x)]s + exp. p(y)dy] = 0 
(17) 


SOLUTION OF THE EQUATION OF FREE OSCILLATIONS 


To construct the general solution of Eq. (6) or Eq. (5), it is therefore necessary to find two particular solutions 


satisfying condition (15). 


Two such solutions are easily obtained. For that purpose 6; and #2 are expanded in a series, progressing in 


powers of the advance ratio yu: 


Bn = Bro + Bart + Baru? +.... (nm = 1, 2) (18) 


This is possible because » < 1. Normally, u will hardly exceed the value 0.3. 
In order to satisfy Eq. (15), the following initial conditions are imposed on the coefficients of the series (18): 


= 1, = 0, 


B2o’(0) = 1, 


= 0, 


Bam(0) = 0 


Bam'(0) = 0 
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Using the abbreviations 


7(B?/8) =a, ka =a, 


1(B*/12) = 


y(B*/16) = (20) 


introducing series (18) in Eq. (5), and equating coefficients of like powers of u, one obtains the following infinite 
system of coupled linear differential equations of the second order in the functions 8, : 


+ 2B am’ + (1 2ke1) Bam (2 — COS 2y sin Ba, m—2— 


2cs sin ¥B'n, m —1 — 


2c3(cos — 2k sin = 1,2; m=0,1,2,3...) (21) 


where By, m—2 = Bam—1=0. 


Eqs. (21) can be resolved successively. Carrying the calculations tom = 2 only—i.e., terminating series (18) 


after the third term—one obtains 


= Axo cos + Byo sin ef 


= cos ey + B, sin ey + E,; cos (1 + + Fu sin (1 + + Gai cos (1 — + Hy sin (1 — ov | 
= cos ef + sin + cos eb + sin eb + Ege cos (1 + + Fre sin (1 + + } (22) 
Gi2 cos (1 — + Hye sin (1 — + cos (2 + + Kyzsin (2 + + Lye cos (2 — ey + 


sin (2 — = 1, 2) 


where, on the basis of Eqs. (19), the coefficients are found to have the following values: 
Aw = By = Axa = 0, Ba => l/e 
[¢s/(1 + 2e) + m—1 + + 2k)Ba, a-t} 
[c3/(1 + 2e) + By, m—1 — (G4 + 
Gam = [¢3/(1 2e) m—1 + 2k) B,, 
Aan = [—«¢3/(1 2e) ©)B,, + (a + 2R)An, 
(m = 1, 2) 
An = + Gm), Bu = —(1/e)(Q ©) Fa + (1 
Cup = {cz + [2cs2/(4e2 — 1)][1 + + 2k))} (23) 
Diz = (Ano/2€) + [2c32/(4e — 1)][1 + + 2k) 
= [1/8(1 + ©) — Bao) + 2cs[(2 + + (ca + 2k) 
Ky = + €) + kByo) + 2¢3[(2 + ©) Fra — (a+ 2k) 
Le = [1/8(1 €) ]fa(kA no + By) + 2c3[(2 — + (cs + 2k)Hm)} 
M2 = {1/81 €) {alAno RByo) 2c3[(2 (cs 
An = —(Exa + Gro + + L,2) 
= —(1/e)[Co2 + (1 + Fae + (1 — + (2 + + (2 — Mie) 
(n = 1, 2) 


CRITERION OF STABILITY 
The characteristic Eq. (17) of Eq. (5) is 
— [Bi(2n) + = 0 (24) 
From Eqs. (22) and Eqs. (23) it follows that 
B10(2m) + = 2e~*** cos 2we (25a) 
Bu (2m) + Bx’(2r) = 0 (25b) 


where 
2 
c= + [1 + + (25d) 
€ de 
Writing 
s = oe (26) 


Eq. (24) transforms into 


o? — 2(cos + wx sin + ...)o +1=0 (27) 


The two roots o; and o: = 1/0; of the characteristic 
Eq. (27) can therefore easily be determined, and there 
is no necessity to determine all the coefficients defined 
by Eqs. (23). 

It follows from Eq. (14) that the disturbed motion is 
stable if 


|s| = <1 (28) 


When ¢ is complex, then |c| = 1. 

In Fig. 1 |s| has been plotted against the advance 
ratio u for different values of the blade mass coefficient 
y when B = 0.98, both for k = 0 (no coupling between 
flapping angle and blade angle) and k = —0.75 (strong 
coupling). It may be seen that |s| remains constant 
for lower values of u. Then, depending on y, |s| be- 
gins to increase from a certain value of u on, indicating 
that partial destabilization sets in. Nevertheless, in 
the u-range shown, |s| always remains far below the 
critical value of |s| = 
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THE SPECIAL CASE, ¢€ = 0 


In this case Eqs. (22) are not valid. Floquet’s 
theory here leads to a solution of the form 
By = Cree“ + (29) 


where yg, and ¢» are periodic functions with the period 


This solution is stable at all u for any given value of 


y because cy > 0. 


EXAMPLE 


A rotor with the data! 


SCIENCES—FEBRUARY,1948 


has been calculated both for » = 0 (hovering condition) 
and for » = 0.3. Noting that’ 


k Ce C3 € 
0 1.44 0 0.93 0.69 0.725 
—0.75 1.44 —1.08 0.93 0.69 1.249 


one obtains the characteristic Eqs. (27) in the form 


k=0: o? + 2(0.156 + 6.708u7)o + 1 = 0 
k = —0.75: o? — 2(0.006 — + 1=0 


The values of |s| = |c\e~?** = 0.013 |o| have been 
plotted in Fig. 1. The flapping motion is stable be- 
cause |s| < 1 for all u considered. 

In determining the general solution, Eqs. (18), it is 
necessary to evaluate all the coefficients, Eqs. (23). 
They are given in Table 1 for k = O (upper line) and 
k = —0.75 (lower line). 

The particular solution 6;, that satisfies the initial 
conditions 6,0) = 1, 6,/(0) = 0, is plotted in Fig. 2 
for one revolution of the blade, both when » = O and 
when » = 0.3. It is seen that the disturbed flapping 
motion is nearly extinguished after one revolution of 
the rotor. 
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A Semigraphical Method of Computing Stick 
Forces for Spring-Tab Controls Having 
Nonlinear Hinge-Moment Characteristics 


ROBERT A. MENDELSOHN? 


Continental Aircratt Corporation 


SUMMARY 


A method of computing stick forces for spring-tab controls is 
presented in which aileron and tab characteristics from wind- 
tunnel tests are used instead of linearized characteristics. To 
illustrate the method, the variation of stick force with rate of 
roll is calculated for two flight speeds for an assumed airplane. 
The calculations are made by the semigraphical method presented 
herein and by a method using linearized aileron and tab data to 
illustrate the desirability of the more complete semigraphical 
method for large rates of roll. 

Using a chosen differential aileron control, a determination of 
the optimum spring strength for high rates of roll with low stick 
force is demonstrated for a high-speed flight condition. Far 
equilibrium, it is shown that the optimum spring strength will be 
dependent on the magnitude of the tab hinge moments, and for 
large and fast airplanes, excessive stick forces may occur for an 
unbalanced tab even if the tab has nospring (servotab). An ex- 
tension of the semigraphical method of solution to the case of a 
geared spring tab is indicated. - 


INTRODUCTION 


— OF THE NECEssITY of keeping stick forces 
within the limits of a pilot’s physical strength, the 
use of spring-tab controls is becoming more and more 
prevalent as airplanes become larger and faster. The 
computation of stick forces for this type of control be- 
comes involved, if linear aileron and tab characteristics 
are not assumed, because of the great number of pos- 
sible variables. A semigraphical method of calculating 
stick forces for spring-tab controls is presented in which 
actual control characteristics from wind-tunnel tests 
may be used. The method permits calculation of nec- 
essary spring strength, takes into account tab hinge 
moments, and accounts for conditions during a man- 
euver. 


SYMBOLS 


The coefficients and symbols used in this report are defined as 
follows: 


CL = lift coefficient (L/qS) 


Ch = hinge-moment coefficient (H/gbc?) 
= rolling-moment coefficient (L’/gSb) 


Received April 14, 1947. 

* The author expresses his thanks to Miss Sadie Miller of the 
Langley Memorial Aeronautical Laboratory for carrying through 
the computations and for her assistance in editing the paper. 

+ Head of Aerodynamics Department. Formerly with 
N.A.C.A. 


Cy = rate of change of rolling-moment coefficient C, with 
helix angle (pb/2V) 

L = lift of complete airplane 

W = weight of complete airplane 

H = hinge moment 

i ig = rolling moment, at plane of symmetry, caused by 
deflection of one aileron 


Fs = stick force 
p = rolling velocity 
pb/2V = helix angle of wing tip 
S = area of wing 
b = span 
c = chord 
c = root-mean-square chord 
l = control-stick length 
q = dynamic pressure 
a = angle of attack 
6 = control deflection relative to surface to which it is 
hinged, positive downward 
0 = control-stick deflection 
a = section lift-curve slope, per deg. ~ 
Aa, = effective change in angle of attack caused by rolling 
velocity 
ky = ratio between angular control-stick deflection and 
aileron deflection, spring tab fixed 
ke = ratio between angular control-stick deflection and 
spring-tab deflection, aileron fixed 
Subscripts 
w = wing 
a = aileron 
t = spring tab 
by = spring tab (plus spring strength) 
Ci = rolling-moment coefficient 
Ch = hinge-moment coefficient 


CHARACTERISTICS OF ASSUMED AIRPLANE 


Characteristics of the assumed airplane, used to illus- 
trate the semigraphical method of computing stick 
forces for spring-tab controls, are as follows: maximum 
speed, 360 m.p.h.; minimum speed, 100 m.p.h.; 
weight, 13,000 Ibs.; aspect ratio, 6.24; taper ratio, 
0.434; wing span, 48.97 ft.; wing area, 383.89 sq.ft.; 
aileron span, 10.31 ft.; aileron root-mean-square 
chord, 1.43 ft.; aileron chord/wing chord, 0.20; 
spring-tab span, 3.17 ft.; spring tab root-mean-square 
chord, 0.43 ft.; maximum stick deflection, + 13 deg.; 
maximum tab deflection (with stops), = 15 deg.; 
maximum aileron deflection with tabs against stops, 
—18, 11 deg.; maximum aileron deflection with tabs 
neutral, —25.7, 18.7 deg.; stick length, 27.75 in. 
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Fig. 1 presents a schematic diagram of the assumed 
aileron and tah linkage and also gives the differential 
characteristics as computed from the Appendix. 


METHOD 


To illustrate the semigraphical method of calcula- 
tion, stick forces have been computed for a high-speed 
and a low-speed flight condition of an assumed air- 
plane. Computations were made assuming aileron 
and tab characteristics, as presented in Figs. 2, 3, and 
4. In addition, computation of stick forces is made by 
assuming linear aileron and tab characteristics as 
measured by the slope of the curves shown in Figs. 2, 
3, and 4.! 

Assumptions 

(1) A differential linkage illustrated by Fig. 1 is 
used. 

(2) Aileron and tab characteristics are given by 
Figs. 2, 3, and 4. 


(3) pb/2V = AC,/(Ci,/a0)a0 (1) 


- [this equation to be used only for determining (Aa,)], 


and 
pb/2V = 0.8 [ACi/(Ci,/a0)a0] (1a) 


where the overall correction to pb/2V for wing twist, 
sideslip, and yawing was taken as 20 per cent, where the 
coefficient of damping in roll C,,/ao was taken as 4.18 
from reference 2, and where dp was assumed to be 0.105 


Assumed linkage for the spring-tab aileron used as an - 


for conditions under which data of Figs. 2, 3, and 4 were 
obtained. 


(5) A value of 0.8 of maximum speed of the assumed 
airplane was used in the calculations for the high-speed 
condition (C, = 0.160; @ = 1.5°), and a value of 1.4 
times the minimum speed of the assumed airplane was 
used for the low-speed example (C, = 0.634; a = 
7.7°). 


(6) . (Aa,)C, = 34.4 (pb/2V) (3) 
(Aay)C, = 26.13 (pb/2V) (4) 


as obtained from reference 1. 

(7) A maximum tab deflection of approximately 
+ 15° will be reached when Fs = 30 lbs. 

(8) The ailerons are assumed to be trimmed near 
a=0°; = 0°. 


Procedure 


_ Because tab hinge moments act effectively as an 
increase in spring strength and because it is desired to 
reach full tab deflection of 15° for the high-speed rolling 
condition to obtain as great an aileron deflection as 
possible for 30 Ibs. of stick force, the determination of 
the required spring strength for the high-speed rolling 
condition will first be attempted. 

The required spring strength for full 15° tab de- 
flection and 30 Ibs. of stick force will not necessarily be 
the same for the upgoing and downgoing wing of a roll- 
ing airplane because of the differential linkage and be- 
cause of the difference in the tab hinge moments. Since 
only one size spring is to be used for both ailerons on the 
airplane, the greater strength requirement will be used 
as a criterion in order to maintain aileron control at low 


speeds. 
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Fic. 2. Variation of aileron characteristics with aileron de- 
flection for various spring-tab deflections. 
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COMPUTING STICK FORCES FOR SPRING-TAB CONTROLS 


For full-stick deflection, and assuming 30 Ibs. total 
stick force, because of the difference in k, between the 
upward-deflected aileron and the downward-deflected 
aileron and also because only one tab may be deflected 
to the full 15°, the contribution of each aileron to the 
total stick force will not be the same. 

As a first approximation, assume that the stick is 
deflected to its maximum allowable value of 13° and 
calculate, by means of Eq. (2), the value of C,, (as- 
sumed equal for upgoing and downgoing ailerons) which 
will cause a total of 30 Ibs. of stick force at high speed. 


(30) (0.160) 
388.2[(1/0.452) + (1/0.874) ] 


For full-tab deflection of 15°, it may be seen from 
Fig. 2 that an aileron deflection of only 16.7° may be 
reached for the upward-deflected aileron and an aileron 
deflection of 11.3° may be reached on the downward- 
deflected aileron. Thus, for the assumption of 30 Ibs. 
of stick force and for the assumption of both tabs de- 
flected 15°, full-stick deflection is limited by the up- 
ward-deflected aileron. Fig. 1 indicates that the stick 


Cra = = 0.00368 


deflection corresponding to 6, = 15°, 6, = —16.7° is 
12.4°. 

Remembering that @ = 12.4° for the critical aileron, 
it may be shown that 6, = 11° and 6, = —13.4° for 


equilibrium of the downward-deflected aileron. Since, 
contrary to the assumptions, full-tab deflection has not 
been reached for the downward-deflected aileron, the 
stick force will be less than 30 Ibs. when the tab on the 
upward-deflected aileron just reaches 15°. An ad- 
justment of allowable stick deflection for 30 Ibs. of stick 
force may therefore be made. A second approxima- 
tion is given by a modification of Eq. (2), which assumes 
that the aileron hinge moment, balanced by a spring 
tab, varies linearly with tab deflection. 


(W/S)baCa? 5; 6: 


This assumption is reasonable and is sufficiently ac- 
curate for the present calculations: 
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Fic. 3. Cross plot showing variation of aileron hinge- 
moment coefficient with spring-tab deflection for various aileron 
deflections. 
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Fic. 4. Variation of spring-tab characteristics with spring-tab 
deflection for various aileron deflections. 


1 + 13.4 | = 
0.455 -15(0.855) 
0.00381 


Cra = (30) (0.160) /388.2 


With this new value of C,,, no large change in stick 
deflection, aileron deflection, or tab deflection occurs for 
a stick force of 30 Ibs. Because the effective angle of 
attack of each wing of an airplane during a rolling 
maneuver is different from the level flight value, the 
previous calculation should be made using data similar 
to Figs. 2 and 3 for several angles of attack to obtain 
the conditions for the rolling wing. Computing in a 
similar manner the aileron deflections and tab deflec- 
tions for several angles of attack, then the rolling- 
moment coefficients are determined from Fig. 2 and the 
rates of roll pb/2V are calculated from Eq. (1). The 
angle-of-attack change, caused by roll (Aa,)C, as de- 
termined by Eq. (3), may then be added to the down- 
going wing and subtracted from the upgoing wing. 
Since the value of (Aa,)C, varies with angle of attack 
because of the change in equilibrium conditions for the 
aileron, the value corresponding to the rolling wing is 
determined by a process of iteration. 

Before a spring having the proper strength may be 
chosen, it will be necessary to determine the effective 
spring strength contributed by the tab hinge moment. 
Assuming the value of effective spring-strength coef- 
ficient (tab hinge moment plus spring-strength effect), 
Cz, to act as a balance for the aileron hinge-moment 
coefficient, account must be taken of the relative size of 
aileron and tab plus the difference in mechanical ad- 
vantage. 


= 1 (6) 
For the assumed aileron and tab: 
Crp = 22.94 Cy, (7) 


Thus, for the stick force caused by tab deflection to 
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Fic. 5. Variation of spring-tab deflection with aileron deflection 
for several angles of attack and for a high-speed condition. 


equal the stick force caused by aileron deflection, 
the tab hinge-moment coefficient C,,; must be 22.94 
times as large as the aileron hinge-moment coeffi- 
cient. 

The tab hinge-moment coefficient of the upward- 
deflected aileron was found from data illustrated by 
Fig. 4 to be —0.020 for the high-speed rolling condition, 
and the required aileron hinge moment was found by 
a modification of Eq. (2) to be 0.00382 for 30 Ibs. of 
stick force. The required spring strength will be: 


388.2 
(0.455) (0.160) E 22.94 


or a spring strength of 1.05 Ibs. stick force per 
degree tab deflection for the high-speed rolling condi- 
tion. 

Having determined the necessary spring strength, 
the aileron deflections and tab deflections corresponding 
to given stick deflections will be determined for the 
chosen flight speed and for various angles of attack. 
This is done by an iteration process as follows: 

On Fig. 3, which is a cross plot of Fig. 2, and for 
6, = 15°, the spring-strength value of equivalent C,, 
is plotted as a dashed line. A linear variation with tab 
deflection will then be assumed, since it will be inde- 
pendent of aerodynamic forces. To this spring- 
strength value of equivalent C,, will be added the 
equivalent C,, contributed by the tab hinge moment 
corresponding to each tab deflection. As a first ap- 
proximation, it will be assumed that the total equiva- 
lent C,,, contributed by the spring and the tab hinge 
moment, varies linearly with 6,. At a chosen tab de- 


flection and angle of attack and for the assumed value 
of C,, from the linear variation of total C,, with 6, a 
corresponding value of 6, is determined from Fig. 2. 
For the value of 6, just determined and for the assumed 
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Fic. 6. Variation of aileron deflection with stick deflection for 
several angles of attack and for a high-speed condition. 


6, and a, a value of Cy, is found from Fig. 4. By means 
of Eq. (7), the equivalent C,,, is determined and added 
to the spring-strength value. If the total value of 
C,, thus approximated does not equal the assumed value 
from the linear variation assumption, a new determina- 
tion of 6, from Fig. 2 may then be made, using the value 
of total equivalent C,, just determined and repeating 
the procedure. By this means, a curve of equivalent 
Cr, against 6, corresponding to a chosen flight condition 
is determined and plotted for each angle of attack, as 
illustrated on Fig. 3. A cross plot of tab deflection 
against aileron deflection from the intersection points 
on this figure and for the high-speed condition is pre- 
sented in Fig. 5 for these angles of attack. Using the 
known aileron deflections and tab deflections from Fig. 
5, the corresponding stick deflections are determined 
from Fig. 1 and are shown plotted in Fig. 6 for each 
angle of attack. 


Curves of stick force Fs against rate of roll pb/2V 


may now be obtained by the use of a data sheet similar 
to Table 1 and by the use of Figs: 4 and 6. It may be 
noted that the data sheet is arranged so that successive 
approximations to stick forces for the rolling wing are 
determined. The method involves determining the 
aileron and tab deflections from Fig. 6 at the flight 
angle of attack. The rate of roll is then estimated from 
the values of C, presented in Fig. 2, and the effective 
change in angle of attack (Aa,)C, is determiried as a 
correction to the flight angle of attack by the method 
of reference 1. The same procedure is continued until 
no change in rate of roll occurs for a given stick deflec- 
tion (see Table 1). Aileron hinge moments for the 


angles of attack corresponding to these rates of roll ° 


may be obtained from Fig. 2, and stick forces may be 
calculated using Eq. (2). : 
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COMPUTING STICK FORCES FOR SPRING-TAB CONTROLS 111 


Since it is difficult to determine hinge-moment coefficients accurately from Fig.:2, stick forces were determined 
from the sum of the forces contributed by the springs and the tab hinge moments for conditions in which the tabs 


had not reached the stops. 


15 


CD) 22.94k; 


RESULTS 


Fig. 7 presents the variation of stick force with rate 
of roll for two assumed flight speeds and for the line- 
arized and semigraphical methods of solution. 

Using the proposed semigraphical method of comput- 
ing stick forces, Fig. 7 shows that the assumption of 
linear aileron and tab characteristics may result in large 
errors in computed stick forces at high rates of roll. 

Calculations indicated that large differences in esti- 
mated spring size will result, depending on whether or 
not tab hinge moments have been taken into account. 
If an unbalanced tab is used, the necessary spring is 
weaker because of the tab hinge moments acting as an 
effective increase in spring strength. Because, for an 
unbalanced tab, a weak spring is indicated at high 
speeds and because, at low speeds, the effective spring 
strength caused by tab hinge moment is low, the best 
spring size for high speeds will not be optimum at low 
speeds. Also, for large and fast airplanes, stick forces 
caused by tab hinge moments may limit the rate of roll, 
even for a tab with no spring (servotab). For best re- 
sults at all speeds, a balanced tab is desired. 


CONCLUDING REMARKS 


In the example used to demonstrate the semigraphi- 
cal method of solution, no mention was made of preload 
or of geared spring tabs. Tabs having a preload in the 
spring present no additional problems and may be 
handled as before, except that the curve of equivalent 
Cyg against tab deflection (Fig. 3) will be displaced on 
each side of the origin. The case of geared tabs, how- 
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Fic. 7. Comparison between ‘the linear and semigraphical 


stick-force results. 


ever, will necessitate a slight change in the method of 
successive approximation used to obtain the curve of 


TABLE 2 
COMPUTATIONS OF STICK FORCE AND RATE OF ROLL FOR VARIOUS STICK DEFLECTIONS 
(HIGH-SPEED CONDITION) 


Py Correction Dus to Effect of Rolling 
3 a 1 R 
4.3 8.4 10.9 1.4 11.6 
2.0 49 8.3 15.0 
(he 4.0094 .0170 0198 
(cas, 20304 +0340 0345 
Cz" -0190 0379 0502 -0556 0563 
1.49 2.97 3.9% 4.35 
8.4 1.8 1.9 
“1.9 3.9 27.2 “11.0 
-9.0 16.4 ~18.0 
& 2.5 6.0 10,0 15.0 
(eis, O17%6 .0218 .0220 +0246 
0086 .0193 .0308 0340 
0178 +0526 0560 -0591 
0405 1275 1345 
Barc 1.39 2.89 4.38 463 
1.06 2.19 3.33 3.51 
Se 8.4 11.1 11.8 1.9 
“1.9 “3.9 “7.3 -11.0 “12.5 
-&, 16.4 -18.0 
2.5 6.0 9.8 15.0 
— 0523 — 
(Aép) cy — 3. — 
Spring .00077 4.00295 
tab 
00026 00013 00057 00140 
( ) 00049 00118 -00193 00283 
% tab Oud 
627 692 836 27% 
(me, 2 6.28 12.61 33.08 
(hs 4.7 10.69 2.% 20.00 25,00 
7.17 13.85 19.04 32.61 58.98 
0324 0673 0953 1021 1077 
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equivalent C,, versus tab deflection. Because, for a 
given flight speed, the tab deflection will depend upon 
the aileron hinge moment, the aileron deflection, the 
spring strength, and the tab hinge moment, the equiva- 
lent C,, caused by spring strength and tab hinge mo- 
ment should be plotted on Fig. 3 along a line of constant 
aileron deflection 6, starting from the tab deflection cor- 
responding to each aileron deflection as determined 
from the gearing. 


Appendix 


CHARACTERISTICS OF ASSUMED AILERON LINKAGE 


Assuming a spring-tab linkage as pictured in Fig. 1 
and a maximum control-stick deflection of +13°, the 
following method was used to obtain curves defining 
the aileron and tab motion as a function of control-stick 
motion: The angular deflections of the control stick 
and push-rod quadrant were assumed to be the same 
and the change of aileron angle with push-rod move- 
ment was assumed to be linear. For full 15° tab de- 
flection, the push-rod, aileron-link deflection would be, 


from Fig. 1: 
(1.0937/2.125)15 = 7.72° 


Designing the linkage such that, for full-tab deflection, 
the maximum upward-aileron deflection will be 18° 
and the maximum downward-aileron deflection will be 
11°, the maximum aileron deflections for tab neutral 


would be: 
25.72° 


(up aileron) 7.72 + 18 = 

(down aileron) 7.72 + 11 = 18.72° 
Letting: 
R = radius of push-rod quadrant 


K 


push-rod movement per unit aileron deflection 
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6 stick deflection 

6 aileron deflection 

6, = tab deflection 

@ = angle of push-rod quadrant link 


with the subscripts: 


0 = zero stick deflection 
“u = upward-deflected aileron 
d = downward-deflected aileron 


R(sin — $0) = — (7.72/15)6;] (1) 

R(sin — sin = K[é, — (7.72/15)6,] (2) 

where ba = + Oand = — 8. By substituting 

the conditions for maximum stick deflection in Eqs. 

(1) and (2) it may be shown that ¢ = 54.12° and 

R/K = 168.6°. The equations of motion are there- 
fore: 

168.6 [sin (54.12 + 6) — 0.811] = 6, — 0.51456, 
168.6 [0.811 — sin (54.12 — 6)] = 6, — 0.51456, 
Assuming a constant tab deflection and differentiating 
with respect to control stick deflection, we obtain the 

mechanical advantage of the aileron as: 


1/kig = 2.942 cos (@ + 54.12) 
1/ki, = 2.942 cos (54.12 — 4) 
and the mechanical advantage of the tab: 
ke = 0.5145 
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Matrix Development of Multhopp’s 
Equations for Spanwise Air-Load 
Distribution 


STANLEY U. BENSCOTER* 


National Advisory Committee for Aeronautics 


SUMMARY 


The Prandtl lifting-line equation is replaced by a system of 
linear algebraic equations that govern a finite set of ordinates to 
the lift distribution curve. Matrix algebra is used for the de- 
velopment of the equations. The matrix of the system of equa- 
tions consists of the sum of two square matrices. The first of 
these matrices is a diagonal matrix containing the physical con- 
stants for the particular wing to be analyzed. The second ma- 
trix contains elements that are standard universal constants ap- 
plicable to all wings. A numerical example is given. 


NOTATION 


A = wing aspect ratio 

ad, = Fourier coefficient 

= wing span length 

= wing chord length 

= effective angle of attack 

= induced angle of attack 

= total wing lift 

= wing lift per unit of length 

= section lift curve slope 

= dynamic pressure (p V?/2) 

= wing plan form area 

= velocity 

induced velocity 

= linear spanwise coordinate 

= angle of attack 

= circulation per unit of length 

= a dimensionless circulation (=21'/bV); (= 1/bq) 
mass density of air , 
= angular spanwise coordinate 

a dimensionless section property (= b/mc) 


FTPR MRE SHS 


INTRODUCTION 


§ bes METHODS OF CALCULATION to be illustrated are 
based u,on Prandtl’s lifting-line equation. No 
attempt is made to develop an analytical solution of 
the equation. The method of solution involves the 
solution of a system of linear algebraic equations. The 
unknowns that are governed by the equations are or- 
dinates to the lift distribution curve. When these 
ordinates are plotted at the appropriate spanwise sta- 
tions, the lift curve may be drawn through the points 
thus established. In order to emphasize the use of 
matrix algebra, the analysis will be limited to a con- 
sideration of the simplest wing air-load case. The 
wing is considered to have no twist of any sort. The 
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lifting line, which may be regarded as the line of aero- 
dynamic centers of the wing sections, is considered to be 
straight and normal to the plane of symmetry of the 
airplane. 


INDUCED VELOCITIES 


Vertical velocities arise at points along the wing 
which may be defined as being “induced” by the trail- 
ing vortex sheet. The development of the formula for 
induced velocities will be stated briefly. A single-line 
vortex extending from to + © is shown in Fig. 
la. If the strength of the vortex is I’, the velocity of 
the air at a distance / is given by the formula! ['/2rh. 
In Fig. 1b a semi-infinite vortex is shown extending 
from the origin to +. The velocity at point a, which 
lies at a distance h from the vortex, is given by the 
formula! ['/4rh. The velocities at points to the right 
and left of point a, which are also at a distance h from 
the vortex, are illustrated in Fig. le. The trailing 
vortex sheet may be replaced by a finite number of 


a 


(a) Infinite Vortex 


T 
41h 
(b) Semi-infinite Vortex 


4 Wa = 
a 


Wa" 
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(c) Induced Velocities for 
Semi-infinite Vortex 
Fic. 1. Line vortices. 
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(a) Linear and Angular 
Coordinates 


(b) Lift Distribution 


Fic. 2. Location of lift ordinates. 


isolated vortices for the purpose of making numerical 
computations and for illustration. 

The section lift / is proportional to the section cir- 
culation I in accord with the Kutta-Joukowsky law.’ 


l= pVT (1) 


Because of this proportionality the spanwise lift dis- 
tribution and the spanwise circulation distribution may 
be used interchangeably in discussion. In the present 


method of calculations the section lift, or spanwise lift 


curve ordinate, will be computed at certain specific 
points. The location of these specific points is shown 
in Fig. 2. The distance y to any point along the wing is 
assumed to be related to an angular coordinate @ by 
the formula, 


y = (b/2) cos 6 (2) 


The spanwise stations 1, 2, ys, etc., are determined by 
using equally spaced values 6;, 62, 63, etc., of the angular 
coordinate. 

The circulation distribution over the semispan is 
indicated in Fig. 3a with four ordinates T, to Ty. In 
the interest of brevity, the computation methods will 
be limited to the determination of four lift curve ordi- 
nates. This leads to 4 by 4 matrices, which are suf- 


ficiently large to illustrate the methods of computation. 
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In Fig. 3b the circulation diagram is divided inté 
four horizontal bands. These are formed by drawing 
horizontal lines through points on the curve at the 
spanwise stations where the I values are located. 
Each band may be regarded as representing the span- 
wise segment of a horseshoe vortex. The height of the 
band is the strength of the isolated vortex. These four 
horseshoe vortices may be drawn with their spanwise 
segments superposed as shown in Fig. 3c, or they may 
be drawn separately as shown in Fig. 3d. From these 
illustrations it is apparent that the strength of a trailing 
vortex is equal to the difference of two adjacent ordi- 
nates to the circulation diagram. If a diagram of the 
strengths of the trailing vortices is plotted, it will ap- 
pear as shown in Fig. 3e. This stepped diagram is an 
approximate representation of the derivative dI'/dé, 
since the variable @ has been considered to vary in 
equal increments. If the ordinates to this diagram 
were multiplied by values of d@/dy, which would be ob- 
tained from Eq. (2), the resulting diagram would repre- 
sent the amount of vorticity per unit width of the trail- 
ing vortex sheet. 

The same concepts may be described in the language 
of continuous functions. In Fig. 4a the circulation 
diagram is shown with an ordinate T at station y and an 
ordinate [ + dI at station y + dy. In Fig. 4b an iso- 
lated vortex is shown at the center of the width dy. 
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Fic. 3. Trailing vortices. 
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As before, the strength of this vortex is equal to the dif- 
ference between the circulation ordinates at either side 
of the width being considered and, hence, is equal to 
df. The values of df at various stations obtained for 
equal space increments dy may be plotted to form the 
diagram shown in Fig. 4c. This diagram represents 
the spanwise distribution of vorticity in the trailing 
vortex sheet. : 

Referring to Fig. 4b, the induced velocity at station 
y’ due to the trailing vortex of strength dI! may be 
computed from the formula for a semi-infinite vortex. 
Let dw indicate the induced velocity at y’ due to the 
vortex 


dr 


— (3) 


dw 

The value of the entire induced velocity w at station y’ 

may be obtained by integrating this equation across the 
wing span. 

w= (1/4r) 

It is well known that this integral may be evaluated by 


assuming the circulation to be given by a Fourier 
series. 


= (bV/2) a, sin (5) 


The coefficient of this series is chosen so that the coef- 
ficients of the series will be dimensionless numbers. 

Substituting Eq. (5) into Eq. (4) and integrating 
gives 


w= (V/4sin 6) >> na, sin nd (6) 


This integration is illustrated in various textbooks.‘ 
Having developed a formula for induced velocities by 
the processes of analytical functions, the remainder of 
the process of solution can be performed by various 
numerical methods. It is convenient to introduce a 
dimensionless circulation y, as suggested by H. Mul- 
thopp.® 
= (6V/2)y (7) 
or 
= Da, sin n0 (8) 


The lift curve ordinate is related to this dimensionless 
circulation by the formula, 
= bay (9) 
where 
q = ‘/2pV? (10) 
Thus, one may obtain a lift distribution for a unit 
value of g by multiplying y values by the span, b. 


T} |U+dr 


(b) 


Fic. 4. Trailing vortex sheet. 


INDUCED ANGLE MATRIX 


A matrix of influence numbers for computing induced 
angles at a finite number of stations from circulation 
values may be developed by using matrix algebra. 
The processes of addition, subtraction, multiplication, 
and division of matrices has been explained in an ele- 
mentary manner elsewhere. As is well known, an 
unsymmetrical case is most conveniently solved by di- 
viding it into symmetrical and antisymmetrical parts 
and solving for each part separately. This is permis- 
sible as long as the lift curve is assumed to be a straight 
line. In Eqs. (5), (6), and (8) the odd terms (m, an odd 
number) contribute to the symmetrical part, while the 
even terms (”, an even number) contribute to the anti- 
symmetrical part. When the two parts are computed 
separately, it is only necessary to determine lift curve 
ordinates over the semispan rather than the entire 
span. Throughout this paper the illustrations will be 
limited to the symmetrical case. 


Eq. (8) may be written four times for the stations corresponding to 4, 62, 43, and 4;. 


1 = ad sin 6, + a3 sin 30, + ds sin 50; + az sin 70; 
72 


a, sin 4. + a3 sin 30. + ds sin 50. + a7 sin 760 


(1la) 
(11b) 
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Ys = a sin 63 + a3 sin 363 + ds sin 503 + a7 sin 763 
v4 = a Sin 0 + a3 sin 30, + as sin 56, + a7 sin 764 


The @ values are 
6, = 22.5°, 0. = 45°, 


(11ce) 
(11d) 


63 = 67.5°, 04 = 90° (12) 


If these values are substituted into Eqs. (11), the circulation values become linear functions of the Fourier coef- 


ficients. 
1 0.38268 0.92388 
0.70707 0.70707 
¥3 0.92388 —0.38268 
1 


In contracted form this matrix equation may be 
written 


[vy] = [S][a] (14) 


In this equation [y] and [a] are column matrices, com- 
monly called column vectors. The matrix [S] is a 
square 4 by 4 matrix. 

The matrix [.S] has certain unusually convenient 
features. Each row may be regarded as a row matrix 
or row vector. The elements of a row vector, or column 
vector, may be called the components of the vector. 
The row vectors of [.S] are orthogonal to each other. 
(If corresponding components of two of the vectors are 
multiplied and the result added, zero will be obtained.) 
The norm of the first three row vectors is 2, while the 
norm of the fourth row vector is 4. (The norm of a 
vector is the sum of the squares of its components.) 
The principal diagonal of a matrix runs downward 
from the upper left-hand corner to the lower right-hand 
corner. The transpose of a matrix is obtained by ro- 


The resulting linear equations may be written as a single matrix equation in expanded form. 


0.92388 0.38268 a 

—0.70707 —0.70707 as (13) 
—0.38268 0.92388 as 
1 —1 a7 


tating the matrix about its principal diagonal and is 
indicated with a prime. If the matrix [S] is postmulti- 
plied by its transpose [.S]’, the result is a diagonal ma- 
trix. 


[SIISY = (15) 


oon 
onoo 


The order of multiplication must not be reversed. 
Since, in general, matrices are noncommutative in mul- 
tiplication, the order of multiplication is always of the 
utmost importance. From this result it becomes ap- 
parent that, with slight modifications, the transpose 
[.S]’ can be converted into the reciprocal, or inverse, ma- 
trix [S]~'. In general, a considerable amount of cal- 
culation work is involved in the determination of the 
reciprocal of a matrix. In the present case the recip- 
rocal matrix may be immediately written down as 


0.38268 0.70707 0.92388 0.5 
= 1] 0.92388 0.70707 —0.38268 —0.5 (16) 
2 | 0.92388 —0.70707 —0.38268 0.5 
0.38268 —0.70707 0.92388 —-0.5 


This matrix was obtained by transposing the matrix 
[S] as given by Eq. (13) and changing the values of 
unity in the last row to values of 0.5. The scalar factor 
1/2 is to be multiplied by every element within the ma- 
trix. The product of a matrix by its reciprocal gives 
the identity matrix. 


[S][S]-? = [S]“[S] = (17) 


Thus, the identity, or unit, matrix plays the same part 
in matrix algebra that unity plays in the algebra of 


ordinary numbers. Multiplying Eq. (14) through by 
[S]—! gives, 


= (18a) 
or 

[S]“ly] = (18b) 
or 

= [a] (18¢) 


This equation gives the Fourier coefficients as linear 
functions of the circulation values. 


The ‘‘induced”’ angle at a wing station is defined as the ratio of the vertical velocity, w, to the horizontal veloc- 
ity, V. Using Eq. (6) the induced angle, 7, may be expressed as 
i = w/V = (1/4sin 0) >> na, sin n0 (19) 
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This is a formula for the magnitude of the induced angle without regard to its sign. The induced angle may be 
expressed at four stations as was done with the circulation values: 


4, = (1/4 sin sin 6; + 3a3 sin 3@, + 5a5 sin 50; + 7a; sin (20a) 
ig = (1/4 sin 62)(a; sin 62 + 3a3 sin 362 + 5as sin 562 + 7a; sin 762) (20b) 
iz = (1/4 sin 63) (a; sin 03 + 3a3 sin 36; + 5a, sin 503; + 7a, sin 763) (20c) 
ig = (1/4 sin &)(ay sin 0 + 3a3 Sin 30, + 5as sin 5, + Zaz sin 76) (20d) 


In order to express these equations as a single matrix equation, it is necessary to define two diagonal matrices [A } 


and [N]}. 


sin 6; 0 0 0 
0 sin 62 0 0 
[4] 0 0 sin & 0 (21) 
0 0 0 sin 04 
1 0 0 0 
0 3 0 0 
0 0 0 7 
Eqs. (20) may be written in expanded matrix form to give 
aq, sin 4; sin 0; sin 30, sin 50; sin 76, a 
12 Sin _1] sin ® sin sin 56 sin 762 3d3 (23) 
iz sin 63 4] sin 05 sin 363 sin 563 sin 76s 5ds “i 
sin 04 sin 04 sin 36, sin 56, sin 70, 
The two column matrices may be written as the product of a diagonal matrix and a column matrix. 
sin 0; sin 0 0 0 
sin 62 0 sin 0, 0 0 
in sin | 0 0 sin 0 i, (24) 
14 sin 04 0 0 0 sin 04 ts 
ay, 1 0 0 0 aq, 
3a3 rool 0 3 0 0 a3 
5ds 0 0 5 0 ds 
7a7 0 0 0 ay 


Using Eqs. (21), (22), (24), and (25), Eq. (23) may be expressed in contracted form as follows: 
= '/s{S][V] [a] (26) 


In order to solve for the column vector [i], it is necessary to multiply through by the reciprocal matrix [A]~'. 
Since [A] is a diagonal matrix, the reciprocal matrix is immediately known. The reciprocal of a diagonal matrix 
is formed by replacing each of its diagonal elements by the reciprocal of those elements. Multiplying through 
Eq. (26) by [A]— gives 
= [a] (27) 
Substitution of Eq. (18c) gives 
= (28) 
This equation expresses the four induced angle values as linear functions of the four circulation values. The four 
coefficient matrices must be multiplied together to obtain a single matrix. This multiplication must be performed 


in three steps, each of which will be shown in order to make the calculations clear. The reciprocal matrix [S]—' 
is first premultiplied by the diagonal matrix [N]._ 


0.19134 0.353535 0.46194 0.25 
1.38582 1.060606 ~0.57402 ~0.75 
[NLS] 2.30970 1.767677 ~0.95670 1.25 (29) 
1.33938 9.474747 3.23358 ~1.75 


This result must now be premultiplied by the matrix [.S]. 
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4 1.46499 0 —0.11208 
~1, 46499 4 —1.68916 0 
[SI 0 ~1.68916 4 —1.57716 
—0.22416 0 ~3. 15432 4 


This matrix has certain properties that serve as a check on the calculations up to this point. Except for the last 
row and column the matrix is symmetrical. The nondiagonal elements of the last row are twice the nondiagonal 
elements of the last column. The diagonal elements are all the same, being an integer that is equal to the number 
of stations on the semispan being computed. Thus, all elements are checked. 


The next step is to premultiply by [A]—'. 


10. 45260 —3.82824 0 —0.29288 
—2.07170 5.65656 — 2.38871 9 
-1 
[Ay 0 1.82833. 4.32956 ~1.70710 
—0.22416 0 —3.15432 4 
The elements of this matrix must be divided by 4 as indicated by Eq. (28). It is convenient to define the matrix 
[M] as 
[M] = (32) 
From Eq. (31), [7] is computed to be 
2.6131 —0.9571 0 —0.0732 
—0.5179 1.4141 —0.5972 0 
0 —0.4571 1.0824 —0.4268 
—0.0560 0 —0.7886 1 
This is Multhopp’s induced angle matrix. It is atma- where 


trix of influence numbers which defines the linear rela- 
tionship between induced angles and circulation values. 
The matrix [1/7] has been computed by Multhopp’ for 
solutions with four, eight, and 16 stations on the semi- 
span for both symmetrical and antisymmetrical cases. 
Munk’ has published the matrix for a symmetrical case 
with five stations and the solution at ten stations for 
both symmetrical and antisymmetrical cases. Because 
of differences in defining the dimensionless circulation 
it is necessary to divide the matrix elements as pub- 
lished by Multhopp by 2, while those of Munk must be 
divided by 4 in order to conform to the definitions of the 
present paper. 

The matrix [M] as defined by Eq. (32) may be intro- 
duced into Eq. (28) to obtain 


= [M][y] (34) 
The vector [i] is a linear matrix function of the vector 
[y]. 
LiIFTING-LINE EQUATION 


Prandtl’s lifting-line equation states that the effec- 
tive angle, e, plus the induced angle, 7, must equal 
the angle of attack, a, at every station along the span. 

e+ti=a (35) 


The effective angle may be expressed in terms of the 
section lift, or circulation values, by the formula 


e = 1/mcq = by/mc = py (36) 


a = b/mc (37) 


The physical section property u will have four values, 
HM; to yy, at the four spanwise stations corresponding to 
6; to 4. The column vector of effective angles may be 
written as 


M Y1 m O O O 
= = 38 
Ms 0 O ¥3 ( 
Ma Y¥4 M4 74 


If the diagonal matrix is represented as [u], the above 
equation may be written in contracted form. 


= [z]ly] (39) 


If Eq. (35) is written at the four stations being con- 
sidered, the system may be written as an equation of 
column vectors. 


1 Qa) 
dy a2 

+ = 40 
3 13 a3 ( ) 
14 a4 


This may be written in contracted form. 
[e] + [¢] = [a] (41) 


The vectors [e] and [7] have been previously expressed 
as linear matrix functions of the vector [y]. Substi- 
tuting Eqs. (34) and (39) into Eq. (41) gives 


+ = [e] (42) 
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or [Bl{y] = [a] (45) 
+ [M]} [7] = [a] (43) 
Let This is a linear system of equations governing the cir- 
[B] = [x] + [7] (44) culation values. The matrix [B] may be written in 
Then expanded form as: : 
(2.6131 + ps) —0.9571 0 —0.0732 
—0.5179 (1.4141 + ys) —0.5972 0 
0 —0.4571 (1.0824 + ys) —0.4268 
—0.0560 0 —0.7886 (1 + us) 
In order to calculate the lift curve ordinates for a given Substituting Eq. (5) gives 
wing, the section properties to must be computed +b/2 
and inserted into the matrix [B]. The matrix [M]isa L = bq S073" X ay sin nO dy (48) 
standard universal matrix, applicable to all wings, and From Eq. (2), 
the section properties of a given wing will contribute dy = —(b/2) sin 6 d0 (49) 


only to the diagonal elements of [B]. The linear equa- 
tions represented by Eq. (45) may be solved for circu- 
lation values, which may be converted to lift curve 
ordinates, for a unit value of g, by multiplying by the 
span 


Tota Lirt Force 


When the circulation distribution along the span of a 
wing is assumed to be given by a Fourier series, as in 
Eq. (5), it is well known that the total lift force is pro- 
portional to the first coefficient a;. This fact is de- 
termined by the following integration: 


Substituting into Eq. (48) and changing the I:mits 
appropriately gives 


L = (—b’q/2) f° sin a, sin nd do (50) 
Performing the integration gives 
L= (3b?q/4)ay (51) 


Eq. (18c) gives the Fourier coefficients as a linear func- 
tion of the circulation values. This matrix equation 
represents a system of scalar equations. The first of 
these scalar equations gives the value of a;._ Using the 
first row of the matrix [S]~', as given by Eq. (16), the 


L= St? ldy = f-$'s? pvr dy (47) formula for a; in a four-point solution becomes: 

a, = '/.[0.38268 0.70707 0.92388 0.5) 
v2 
$ ¥3 (52) 

The lift force becomes: 

L = (rb’q/8) [0.38268 0.70707 0.92388 0.5] [x] (53) 


A similar formula for a solution using any number of stations is easily develo ed. 


NUMERICAL EXAMPLE 


A wing of tapered plan form, as shown in Fig. 5, 
will be considered. The lift curve slope will be as- 
sumed to have a linear variation from 5.6 at the root 
to 5.7 at the tip. Linear formulas for the chord and 
lift curve slope may be written as: 


c = 10 — (5y/25) (54a) 
= 10 — 5cos 0 (54b) 

m = 5.6 + (0.1y/25) (55a) 
= 5.6 + 0.1 cos 6 ” (55b) 


From these formulas the four values of the section 
property » may be computed as shown in Table 1. 


TABLE 1 
6 cos 6 m mc 
1 22.5° 0.92388 5.38060 5.6924 30.6285 1.6325 
0.70707 6.46465 5.6707 36.6591 1.3639 
3 67.5° 0.38268 8.08660 5.6383 45.5947 - 1.0966 
5 0.8929 


6 56 


The values of u, as given in Table 1, may be substi- 
tuted into Eq. (46) to obtain the coefficient matrix 
[B]. Assuming the wing to have a unit angle of at- 
tack, the equations appear in expanded matrix form as 
follows: 


30) 

last 
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4.2456 —0.9571 0 —0.0732 a 1 
—0.5179 2.7780 —0.5972 0 (56 
0 —0.4571 2.1790 —0.4268 56) 
—0.0560 0 —0.7886 1.8929 v4 1 
L = 0.6618ASq (59) 
io’ fs : Dividing by Sq gives the lift coefficient as 
C, = 0.6618A = 4.412 (60) 
| 25 | 
CONCLUSIONS 


(a) Planform of Semi-span 


5.6 m 5.7 


(b) Slope of Lift Curve 


0 0.5 1.0 
b/2 


(c) Circulation Distribution 
Fic. 5. Tapered wing. 


The solution of these equations gives the following di- 
mensionless circulation values: 


v1 = 0.394, ¥2 = 0.593 (57) 
= 0.750, = 0.852 
These values are plotted in Fig. 5. 


If the elements of the vector [y], as given above, are 
substituted into Eq. (53), the total lift force on the 


wing may be computed. The product of the row vector 


by the column vector is obtained from a single continu- 
ous machine operation. 


L = (xb’g/8)(1.685152) = 0.66186 (58) 


In terms of area and aspect ratio, this lift may be ex- 
pressed as 


In the language of matrix algebra, a simple procedure 
has been presented for establishing a system of linear 
algebraic equations that govern the spanwise lift dis- 
tribution for an airplane wing. These equations define 
a finite set of ordinates to the lift distribution curve in 
terms of a finite set of ordinates to the angle-of-attack 
distribution curve. The equations were originally de- 
veloped by H. Multhopp. The equations provide the 
natural direct approach to the determination of lift 
distribution for structural design purposes. . 
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ABSTRACT 


Significant improvements in time, cost, and operating tech- 
nique for solving flutter determinants through the use of punched 
card machines are explained. These advances are due to a new 
approach to the handling of complex numbers, to better operat- 
ing techniques, and to superior equipment now available. 


INTRODUCTION 


he computing problems connected with evaluating 

high-order determinants, particularly with regard 
to the flutter problem, are too well known to require 
elaboration. Leppert! explains a method developed by 
the author’for utilizing automatic punch card equip- 
ment for evaluating the flutter determinant, based on a 
trial-and-error method whose mathematical basis is ex- 
plained by Bleakney.* * 


COMPARISON OF METHODS 


The method of reference 1 has these shortcomings: 

(a) It requires a large amount of I.B.M. equipment, 
since the following machines were used in the process: 
(1) two type 601 multiplying punches (one multiplying 
and the other proving); (2) one type 405 alphabetic 
tabulator; (3) one type 517 summary punch; (4) one 
type 513 reproducing punch; and (5) one type 080 
sorter. * 

(b) The true vector nature of the problem was de- 
stroyed by the device used to eliminate division. This 
means that the numbers worked with in the “starring” 
process did not make physical sense. 

(c) Because division was eliminated, the size of the 
numbers changed considerably during the “starring” 
process. This made it necessary to evaluate manually 
each determinant in each order of reduction to deter- 
mine whether the numbers should be offset to accom- 
plish multiplication or division by 10 in order to keep the 
numbers within range. 

The method to be described has none of these short- 
comings and has additional advantages: 

(a) The following I.B.M. equipment is used in the 
process: (1) one type 602 computing punch; (2) one 
type 513 reproducer; and (3) one type 080 sorter. 

(b) The numbers are handled as true vectors. 

(c) The size of the numbers does not have to be con- 
sidered at any point. 
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(d) Complex numbers, as such, are eliminated. 
This means that there are no special sign considerations 
or cross-multiplication and summation problems. In 
addition, although division of complex numbers is ac- 
complished, there is no need for multiplying through by 
the conjugate of the divisor. 

(e) The process is self-checking. 

(f) The operating procedure is so simple that com- 
plete instructions for the reduction could be written on 
a single tabulating card. 


DISCUSSION 


The mathematical basis is exactly as described in 
references (1), (2), and (3). However, the computing 
process is completely different in appearance. If the 
determinant 


= 


whose typical element is a complex number of the form 
a + bi, is rewritten as a determinant of twice the order 
with only a real number in each spot, 

| | 


—bx 


= 


its evaluation will give the determinant squared. A 
slight modification of this principle is used for the flutter 
problem. Instead of placing all of the real elements to- 
gether and all of the imaginary numbers together, each 


- element of the original determinant is split into a two- 


by-two array. For example, suppose that in the orig- 
inal matrix the a, and a): spots were: 


71 72 
| er fe 
27 Zi 


In the double order matrix these would appear as: 
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Working now with the double-order matrix and oper- 
ating as if all terms were real, the ordinary “starring” 
process makes it possible to evaluate the determinant 
correctly. The process is completed when reduced to 
the second order. 

Several items are worthy of note here. First is the 
fact that, by increasing the card volume and eliminating 
operations, the problem has been placed in a form better 
suited to the I.B.M. equipment, which is primarily de- 
signed to handle volume. Secondly, after every two 
reductions—i.e., whenever the order of the matrix is 
even, in each four-element array, representing a single 
complex number in the original matrix—the diagonally 
opposite terms are equal. Since the calculations were 
made in separate cards under different conditions, the 
proof is as complete and elegant as is mechanically and 
mathematically possible. By running the cards 
through a type 070 collator, this comparison can be 
made at the rate of 240 cards a minute. If any errors 
exist, the machine will automatically select these for 
audit. 

The utilization of an International Business Ma- 
chines type 602 computing punch is the second major 
factor in simplifying the flutter problem. This is a new 
machine, which has only recently become available. 

It is a high-speed, relay, counter type of computor, 
capable of the four basic operations (addition, subtrac- 
tion, multiplication, and direct division) and has an 
internal storage and memory device that makes possible 
an entire sequence of calculations. The potentialities 
of this device in the engineering field are undoubtedly 
going to lead to many important developments in com- 
puting techniques. While slow by electronic standards, 
it is available and relatively inexpensive. 

As an example, in constructing the flutter determi- 
nant from the terms 


PA — m+ mQz + mOQgxi 


the machine is capable of remembering the multiplying 
constants P, z, and x and reading from different cards 
the value of the aerodynamic matrix (A), mechanical 
matrix (m), stiffness matrix (mQ), and damping matrix 
(mQgi); automatically determining what matrix terms 


figure. 


are to be multiplied by what constants; summing the 
products; and punching the final determinant amount 
for each element into a trailing.card. To set up a 
particular trial solution of an eighth-order complex 
determinant, engineering specifies a value of P and z, 
a solution order, and the basic A, m, mQ, and mQg 
matrices to be combined. The matrices are identified, 
or called out, by a simple digit code. 

Utilizing the 602 computing punch, the necessary 
matrices can be isolated, combined, and multiplied, and 
the final results for the double order (sixteenth order for 
an original eighth order complex determinant) can 
automatically be punched into the 256 cards required to 
start the reduction process in about 10 min. Trans- 
forming the order of the matrix, to put it in the best 
form for solution, requires another 10 min. The com- 
plete reduction of the determinant by the “‘starring”’ 
process requires about 2'/, hours. If several problems 
are being handled simultaneously, this reduction time 
can be cut to 1'/, hours without any additional man 
power or equipment. 


The actual time for reduction of the determinant is | 


about one-half that required by the method of reference 
1, and the cost is less than one-third of the previous 
If the entire problem is considered (construct- 
ing the determinant, evaluating it, and computing 
amplitude ratios), these time and cost figures become 
even more impressive. 

It should be noted that the system of using a double- 
order matrix can be handled on a type 601 punch, as 
well as on the newer type 602. While a little more 
difficult and time-consuming, it is still superior to the 
previous method. 
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A Simplified Method of Obtaining Drag of 
a High- Speed Body from Wake Surveys 


CHIEH-CHIEN CHANG? 
The Glenn I. Martin Company 


SUMMARY 


Jones’ drag integral equation is re-examined in the light of the 
known experimental fact that nearly constant stagnation tem- 
perature exists across the wake of a body moving at high speed 
in air. It is found that the integrand of the equation can be 
separated into the product of two functions—namely, one func- 
tion containing the ratio of the static to total head pressure at the 
measuring station in the body’s wake and the second function 
containing two similar pressure ratios that correspond to condi- 
tions far up and far downstream along the same stream line. 
Each of the two functions can be conveniently plotted in a graph. 
With the aids of such graphs, this analysis makes the pitot- 
traverse method practical and easy to apply in determining drag 
of flight at subsonic speed. An example of application is shown: 
Some justifications of the assumption are made, and the limita- 
tions of application are explained. 

It is further found that the method can be applied approxi- 
mately to flow around a two-dimensional body with shock and 
possibly to supersonic flow, as long as at and behind the measure- 
ment station no shock waves of appreciable strength exist. Of 
course, some development work has to be done before this 
method can be successfully applied to the supersonic case. 


INTRODUCTION 


i gee PITOT-TRAVERSE METHOD has recently become 
increasingly important in drag measurement of 
flight at subsonic speed. The Jones formula! derived 
from incompressible flow becomes less accurate in such 
high-speed tests where the air compressibility plays a 
significant part. Silverstein and Katzoff? and Lock, 
Hilton, and Goldstein’ revised the formula to take care 
of*the effect of air compressibility. Their resulting 
equations and computation procedures are not quite 
convenient for routine analysis of such tested data. 
These facts and the need of some reasonably simple 
method in a flight-test program have encouraged the 
author to attempt the present analysis. 


THE PROBLEM AND ASSUMPTIONS 


The problem is to find a simple analysis for calculat- 
ing the drag of a body at high flight speed by means of 
pitot-traverse measurement in wake. The assumptions 
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are: (1) Stagnation temperature is constant along a 
stream line; (2) air behaves as a perfect gas; (3) 
the velocity fluctuations in the wake are small in com- 
parison to the local mean velocity that can be calcu- 
lated from the measured time mean pressure; and (4) 
the flow in the wake is isentropic along the stream line, 
downstream of the measuring station. 


DRAG PROBLEM 


From the known laws and relations in one-dimen- 
sional gas dynamics under steady condition, we intro- 
duce here a few equations that will be used in the analy- 
sis: 


dy = (yRTy)”* (1a) 


2 1/2 


b/H = [1 (vy — (Ic) 


} 


2 


y¥-1 


Typ 
(14 


HT (le) 
Here a and dy are the local and stagnation sound veloc- 
ity, respectively; p and H are the static and stagna- 
tion pressure, respectively; Jy is the stagnation tem- 
perature; p is the mass density; g is velocity; / is 
Mach Number; and y is the isentropic exponent, 
which is equal to 1.4 for air. All are referred to the 
same point and based on the isentropic condition. 

For a two-dimensional body in a uniform air stream 
at high speed, as shown in Fig. 1, drag exists, which ac- 
cording to the Jones equation,' can be expressed as 


= — q)ds2 (2a) 


where the subscripts 0 and 2 are referred to the free 
stream and far downstream where, by definition, p, = 
po and ds, is the width of infinitesimal stream-line tube. 
Of course, g is parallel to go. This equation is derived 
from Newton’s second law under the condition of 
steady flow. In the wake, the actual flow is not steady 


123 


the 
__| 
ya 
1Qg 
ed, 
ary 
ind 
for 
can 
l to 
ns- 
est 
ym- 
ig” 
ms 
me 
lan a 
t is 
nce 
ous 
ing 
as 
ore 4 
the 
the . 
‘ical 
ysis, 
-63, 
ctor 
ices, 
| 
= 


124 


JOURNAL OF THE AERONAUTICAL SCIENCES—FEBRUARY, 1948 


SHOCK WAVE 


FAR AHEAD 


Po = const 
Hg = const 
const 


Go = const 


Fic. 1. 


but quasi-steady. The above equation can be shown 
to be true! if we take time-mean value. Thus, 


d = Sor — (2b) 


It can be further simplified by applying the law of 
conservation of mass along the same time-mean stream- 


line tube, pogodso = = pigidsi; thus 
d= — (2c) 


where the station 1 is the measurement station and is 
conveniently located in the wake. It is interesting to 
note that in Eq. (2c) the direction of the measuring 
rack of total head tubes should be parallel to the general 
direction of g, so that dy, can represent the width ds, 
of the stream-line tube under consideration. Until 
now, there has been no such instrument that can meas- 
ure g and pq directly. It is a common practice to 
measure the above quantities by means of pitot tubes 
and to calculate them with the equations given in Eqs. 
(la-le). Because of the time lag of the pressure meas- 
uring system in use, the pressure measured is usually 
the time-mean value. -Therefore, the time-mean values 
of g and pq can be satisfactorily calculated if the velocity 
fluctuation is small in comparison with local mean 
stream velocity. 

Introducing the nondimensional sectional drag coef- 
ficient C, and dropping the average signs, the integral 
can be written as 


(p0/2)qo°C pogo” \qi Jo 


(2d) 


In practice, the integrand is evaluated at each point 
of the wake, and C, is determined graphically with a 
planimeter from the result curve. Thus, we define the 
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TEST STATION 


2 
FAR DOWNSTREAM 


Pi + Po 
: P2 © Po assumed 
Hy + const Hp = Fy along the seme 
tr 


Diagrammatic sketch of the flow pattern around a two-dimensional body at high subsonic speed. 


integrand C,’ as 
dC. (=) ( ( #) 
C,’ = = 2) ——} 
d(yi/ c) Pogo” qo 

2 Bi Qy 

ay ay 
where dy = (yRT x)’ * is constant along a stream line 


under the assumption 1. With Eqs. (1d) and (le) and 
the assumed relation H, = [H,, we may write 


po 
| (=) 
(y¥-D/ (y¥-1)/y7'/2 
)-b-@) 


This expression is a function of the three pressure ratios 
po/Ho, and po/Hi(= p»/H2), which can be ob- 
tained from the wake survey data and the pitot tube 
used for measuring flight speeds. All such data should, 
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line 
and 


Fic. 2. fi vs. 


of course, be corrected for errors inherent in the measur- 
ing system. ; 

Previous papers***® have plotted or tabulated the 
above function as a three-parameter family, which 
makes interpolation tedious. However, examination of 
the above function shows that it can be treated as fol- 
lows: 


Cd’ = «foe (po/Ho, po/ (3b) 


where 
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Fic. 3. foz vs. Ho/po at various 


Here, the integrand has been split into the product of 
two functions, one a single-parameter family and the 
other a two-parameter family. 

Fig. 2 shows f; as a function of p,/, ranging from 0 
tol. Fig. 3 shows fo: as a function of Ho/p) and H,/po, 
each ranging from 1 to 2. Because of the measurement 
setup of the drag problem C,’ cannot be negative. Of 
course, the same equation can be applied to measure the 
thrust of a propeller, where fi. will be negative. 


Instead of po/Hy and po/Hi, Ho/po and po are 
used because greater engineering accuracy ¢an be ob- 
tained for computation at small values of p)/H) and 
po/H;. Fig. 4 shows an actual application of this 
method to flight-test data. From the practical point 
of view, this method is much simpler in derivation and 
application than those of its predecessors. It has been 
extensively used in flight-test work at The Glenn L. 
Martin Company. 
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Drag distribution of NACA 0012 airfoil from wake measurement in flight test. 
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With the transformation relations in Eqs. (la-le), 
Ci’ can be expressed as a function of the Mach Num- 
bers at three stations, but the resulting equation is not 
simple, as can be verified by the reader. However, 
sometimes it is necessary to calculate C,’ for the pres- 
sure ratios beyond the range in Fig. 3, and it is of a 
simple form to express C;,’ in terms of both Mach Num- 


ber ratio and pressure ratios—namely, 


where the Mach Number can be graphically ex- 
pressed asa function of p/H, as shown in Fig. 5. 

In practice, it has been found that the sensitive altim- 
eters are well suited for measuring values ‘of the ab- 
solute pressures. In such cases, p/H can be graphically 
expressed as a function of the pressure altitudes of the 
total head and static reading of the pitot tube. 


PossIBLE APPLICATION TO SUPERSONIC FLOW 


If the above method is applied to the supersonic 
case, the total head readings are no longer correct if the 
local Mach Number is larger than 1, because shock wave 
necessarily exists ahead of the total head opening. For 
the case that the total head tube is in the direction 
local stream line, Lord Rayleigh’ first gave the relation 
between the two pressure ratios p/H’ and p/H, where 
H’ is the measured total head reading and H is the true 
one. Thus, we have 


{ + 1)/2] M?} 


(4 
++ 1 +.1 
+ 
(2 — x 


— 1) 


and briefly 


H’ 
(7M? — 


46,656[(H/p)* — 1] 

{35 — - 

when M > 1 or when H’/p > 1.893. Of course, when 

M<1,H’ =H. 

At the present time, no experimental data are avail- 
able about the Eq. (6a) in the case of the total head tube 
making a small angle of local stream line. Therefore, 
it is advisable to keep the total head tube as parallel 
to the general stream direction as possible. Further- 
more, whether Eq. (6a) can be applied to the wake with 
local large velocity gradient or not, some experimental 
verification is necessary. Some future research work is 
necessary before this method can be applied success- 
fully in the supersonic case. 


DISCUSSION 


The first assumption can be reached by considering 
the air unable to transfer heat across stream line. But 
this is not the actual phenomenon. The assumption is 
approximately verified by the experiments carried out 
by Hilton.* Theoretically, in the wake, two important 
regions may be approximately divided: namely, the 
region that is influenced by the boundary-layer wake 
and the region that is influenced essentially by the 
shock wave, as shown in Fig. 1. In the boundary- 
layer wake, the behavior can be traced back to the 
properties of the boundary layer. In the case of laminar 
boundary layer with zero pressure gradient, Busemann™ 
proves that the stagnation temperature is constant 
across the boundary layer under the assumption of 
Prandtl Number being unity. This simply means that 
inside the boundary layer the rate of production and the 
rate of transmission of heat just balance each other, 
because the ratio of the two rates is the Prandtl Num- 
ber. Actually, the Prandtl. Number of air is 0.733 
under standard conditions. Later, Emmons and 
Brainerd’ show that, for Mach Number = 1, the stag- 
nation temperature is nearly constant across the lam- 
inar boundary layer for actual Prandtl Number of 
air being 0.733. In most tests on airfoils, the boundary- 
layer wake is most probably turbulent. Squire* shows 
that the stagnation temperature across a turbulent 
boundary layer with zero pressure gradient is also con- 
stant under the same assumption of Busemann. For- 
tunately, according to late information, in the turbulent 
boundary layer, the Prandt] Number of air is approxi- 
mately equal to 1, except the laminar sublayer where 
the Prandtl Number is still 0.733. Therefore, the 
theoretical development is substantially conformal to 
the experimental results as far as the stagnation tem- 
perature in boundary-layer wake is concerned. 

In the wake behind the shock wave in the transoni¢ 
case, the viscosity does not affect the flow because the 
velocity gradient there is so small. Besides, the tem- 
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perature gradient is also small, and, consequently, the 
heat transfer is negligibly small. The energy equation 
(q°/2) + GT = constant = C,Ty must be true inside 
that zone of the wake as well as outside the wake. 

In this range of Mach Number, assumption (2) is 
nearly true as far as the specific heat coefficients and 
isentropic constants are concerned. But the mois- 
ture content of the air may cause condensation shocks, 
which is an uncertain implication to the determination 
of the drag. Therefore, the test must be carried out in 
air that has low humidity so that the condensation 
shock surely will not occur. As far as the assumption 
(3) is concerned, the pressure measurements cannot 
give instantaneous readings but give time mean read- 
In the wake, the flow is not steady but quasi- 
steady. Unless the velocity fluctuation is small in 
comparison with the free-stream velocity, the measured 
mean pressure is always slightly higher than that corre- 
sponding to the mean velocity.. 

In cases of separation in the wake, the velocity fluc- 
tuation may become the same order of the mean veloc- 
ity; the measured pressure will be far higher than the 
time-mean pressure. Therefore, the pitot-traverse 
method cannot be applied to the case of flow separation 
without doubt in accuracy. 

Assumption (4), H; = HM, is entirely hypothetical. 
If we can measure q at far downstream, a part of energy 
will no longer be maintained as kinetic energy but will 
dissipate into heat that cannot be measured by means 
of pitot tube. As pointed out in references 6 and 12, 
the Jones formula is more nearly correct only when the 
measurement is carried out near to the body tip where 
the heat dissipation is negligible. 
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Services 


of the Libraries of the 


Institute of the Aeronautical Sciences 


The services of the Libraries are available to all mem- 
bers of the Institute, to Corporate Members, to advertisers 
in the ENGINEERING Revizw and 
NAUTICAL ENGINEERING CaTALoG, and, under usual library 
limitations, to the public. Four specialized services are 
available. 


The Paul Kollsman Lending Library 


This lending library service makes available, without 
charge, the latest and more important aeronautical books. 

Members may request the loan of any aeronautical or 
technical book they wish to borrow. Through an exchange 
agreement with the Engineering Socictics Library, any 
book on general engineering may be borrowed from its great 
collection of over 160,000 volumes. 

A photostating service is available at usual library rates. 

Applications for membership in the library and further 
information will be sent on request. 


The W. A. M. Burden Reference Library 


This reference library contains over 12,000 aeronautical 
books, magazines, pamphlets, and reports gathered from 
world-wide sources and is one of the most complete aero- 
nautical libraries in the world. Material from this library 
is not available for loan but may be used for reference pur- 


poses. 


The Pacific Aeronautical Library 


6715 Hollywood Boulevard 
Los Angeles 28, California 


Established in cooperation with the aircraft companies 
the library serves. The leading aircraft companies in or 
near Los Angeles participate in its support and operation. 

This service library for aeronautical research is available 
to the public for reading privileges. Source material in- 
cludes aerodynamic and structural research reports, as well 


as books on drafting, production methods, history, and al- 
lied sciences. It furnishes books, periodicals, and pamphlet 
material to the participating aircraft companies to supple- 
ment their engineering libraries. 


Technical Information Service. 


This service has experienced personnel under the super- 
vision of trained aeronautical engineers to compile any in- 
formation desired. The services range from listing special- 
ized reference books to the preparation of exhaustive 


bibliographies, digesting of reports, and general surveys of | 
any acronautical subject. Some of the available services | 


are: 


Bibliographies on any aeronautical subject. 

Reports on any aeronautical] subject. 

Digests of aeronautical books, papers, periodicals, and refer- 
ences. 

Translations. 

Engineering investigations of special aeronautical subjects. 

Biographies of individuals engaged in aeronautics. 

Photostats of any aeronautical or general enginecring 
material. 

Microfilms made on special order. 

Photographs made from the Institute’s photographic collec- 
tion. 

Drawings and tracings made. 


In addition to the services mentioned any commission 
which comes within the scope of the Service will be ac- 
cepted. Special arrangements may be made for work* re- 
quiring several weeks or months. 

Translators are available for accurate transcriptions of all 
foreign language data. Translations are carefully edited 
by trained engineers. 

Reproductions of any material in the collections of 
the Institute may be ordered at standard photostat rates. 
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